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POSSIBLE EXPLANATION OF THE “COEXISTENCE” OF FERROMAGNETISM 
AND SUPERCONDUCTIVITY 


B. T. Matthias and H. Suhl 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received December 17, 1959) 


A number of dilute solid solutions of rare earths 
(particularly gadolinium) in host superconductors 
are apparently capable of superconductivity and 
ferromagnetism at the same time.'»? A typical 
diagram of superconducting and ferromagnetic 
transition temperatures, T,; and T,, versus rare 
earth concentration is shown in Fig. 1. For low 
concentration T; >T,, and for higher concen- 
tration T,, > T;. In the shaded region the super - 
conductivity and ferromagnetism appear to co- 
exist. This coexistence is at first sight hard to 
reconcile with the BCS theory of superconduc- 
tivity, * even though the same theory satisfactorily 
accounts for the initial decline of T, versus con- 
centration at low concentrations.* Also there 
occur a number of puzzling experimental effects. 
The object of this note is to propose a simple 
explanation which appears to account for the 
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FIG. 1. Superconducting and ferromagnetic tran- 
sition temperatures, T, ‘g and To as functions of rare 
earth concentration. 


experimental mysteries, and is also consistent 
with the BCS theory. The suggestion is that the 
superconducting regions extend only through the 
thicknesses of the ferromagnetic domain walls. 
Such an explanation is consistent with the BCS 
theory for the following reason. The ferromag- 
netism is almost certainly the result of indirect 
rare earth spin interactions via s -f exchange 
with conduction electrons. A spatially uniform 
spin alignment is difficult to reconcile with BCS 
theory, because it would demand (by virtue of 
the same s-f exchange mechanism) a net uniform 
polarization of the conduction electrons. It is 
possible for such polarization to occur in the 
BCS state because the electron energy in the 
average exchange field exceeds the gap energy. 
However, the normal state permits much greater 
lowering of the energy by such uniform polari- 
zation, and in fact becomes energetically much 
more favorable. On the other hand, in a domain 
wall, where the direction of the rare earth spins 
changes rapidly within one coherence length, the 
effects of the s -f exchange on the energies of 
normal and superconducting states are practi- 
cally equal and so the superconducting state 
remains lower. This is a consequence of the fact 
that the electron susceptibility, in response to 
a spatially varying exchange field, becomes 
obviously the same for normal and superconduct- 
ing states, when the variation is rapid within one 
coherence length. 

We now quote some of the experimental evi- 
dence to support the above viewpoint. Consider 
first the case where the concentration is in the 
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range in which T,>T,;>0. The sample is mag- 
netized at a temperature T, where T,>T>T;, 
and the magnetizing field is removed. The re- 
manence is noted. The temperature is then 
reduced to T<T,;. It is now found that the dif- 
ferential susceptibility equals that of a super- 
conductor. Yet the remanence (as measured by 
removing the sample from a coil and noting 
the galvanometer deflection) is equal to the 
remanence in the previous case T, >T>T ;. In 
terms of the above suggestion, in the range 
T<T,, superconducting regions, coincident 
with the domain walls, form an intricate honey- 
comb or “sponge,” which will obviously give 
the same differential susceptibility as a bulk 
sample. On the other hand, the remanence is 
due to those magnetization vectors M that have 
a component normal to the surface. At those 
surface elements where M does have a normal 
component there is, in general, no domain wall, 
no superconductivity, and thus no difficulty in 
observing essentially the same remanence as in 
the range T,,>T> Ts: 

The full magnetization curve measured by 
Bozorth® on Gd). ogC€o.g2.Ru, (which satisfies 
T, >T >T,) is consistent with the above view. 
Even at the highest fields (~ 12 kilo-oersteds) 
a differential susceptibility characteristic of the 





superconductor still remains. This may be 
accounted for by the remaining ferromagnetic 
domain walls. If these are also superconducting 
and still moderately numerous, a superconduct- 
ing differential susceptibility should remain. 

Next consider the case where the concentration 
is in the range T, >7,,>0. Here we have some 
preliminary results of Phillips® which for Gd in 
La suggest a spin ordering at 7.<T,;. He ob- 
serves a specific heat peak consistent with ferro- 
magnetism. Our picture suggests that the ma- 
terial, fully superconducting for Ts >T> To» 
turns into the sponge-like configuration of non- 
superconducting ferromagnetic domains, sepa- 
rated by superconducting domain walls. 





‘Matthias, Suhl, and Corenzwit, Phys. Rev. Letters 
1, 449 (1958). 
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Solids (to be published). 
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5R. M. Bozorth, Proceedings of the 1959 Conference 
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OBSERVATION OF THE ENERGY GAP BY LOW-TEMPERATURE 
PENETRATION DEPTH MEASUREMENTS IN LEAD 


M. P. Sarachik, R. L. Garwin, and E. Erlbach 
International Business Machines Watson Laboratory, Columbia University, New York, New York 
(Received October 30, 1959; revised manuscript received December 23, 1959) 


The existence of a gap in the energy spectrum 
of electrons in a superconductor implies that all 
observables must become constant at 0°K faster 
than any power of the temperature. Thus, 
exp(-a TJ T) behavior has been observed in the 
electronic specific heat,* in the nuclear spin 
relaxation in superconductors, ? and in ultrasonic 
attenuation, * but not thus far in the contribution 
of “normal electrons” to the electrical losses, 
nor in the penetration depth. By a sensitive 
measurement of the penetration of 2.2-Mc/sec 
magnetic fields through thin films, we have been 
able to demonstrate the effect of the energy gap 
on the low-temperature behavior of the penetration 
depth of lead, a metal in which exponential be- 
havior of the specific heat has not yet been ob- 
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served.* Our data show an energy gap greater 
than 4.9kT, which is considerably higher than 
the 3.527, given by simple theory. 

Our method? is similar to that of Schawlow® and 
Jaggi and Sommerhalder,’ but is considerably 
more sensitive. In Fig. 1, the superconductor 
is seen in the form of a tubular film of thickness 
less than the penetration depth. It serves as a 
shield between the niobium transmitting coil and 
the copper receiving coil, both of which are 
wound in “space harmonic” fashion to produce 
a magnetic field which falls off rapidly with axial 
distance from the coil. This arrangement keeps 
screening currents away from the ends of the 
film and allows shielding factors > 10*® between 
film-out and film-in conditions. The mutual 
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FIG. 1. Diametral section of sample and coils. The 
receiver and transmitter coils are both wound with 
alternate turns in opposite senses to produce and de- 
tect a space harmonic field. The solenoid is used for 
magnetic field investigations which are not discussed 
in the present paper. 


inductance between the two coils is made part of 
an rf bridge, the primary balancing element of 
which is an accurately constructed wave-guide - 
beyond-cutoff piston attenuator (capacitance 10-** 
to 10~'* farad). The bridge is energized by a 
crystal-controlled 2.2-Mc/sec signal, square- 
wave modulated at 400 cps. A 2.2-Mc/sec un- 
modulated homodyning signal and a 400-cps 
phase-sensitive detector provide linear null 
detection with a noise bandwidth of 0.1 cps at the 
signal frequency. In addition to high sensitivity, 
the apparatus has very good stability, since the 
bridge is composed entirely of passive elements. 
With this arrangement, changes in transmission 
of 0.05% can be measured. 

The lead films are evaporated from a multiple 
vapor source onto rotating glass tubes filled with 
liquid nitrogen. The amount of lead evaporated 
for each film determines approximately the thick- 


ness of the films, which is also measured by an 
x-ray fluorescence method.® 

The measurements of the penetration of the 
field through a thin film of thickness d can be 
related to the penetration depth as follows. If 
the correlation length ~ is much less than the 
penetration depth, then the transmission is pro- 
portional to the square of the London penetration 
depth, 27 (a, is the penetration depth in the limit 
&~0.) This result is valid for the London theory, 
the Pippard®,*° nonlocal theory, and the BCS"! 
theory provided that d<< ,;. If d is more nearly 
equal to X;, then one can show that 

2 2 
1 AL (T) es es (0) 
A(T) X pO r 0) ( *) ‘ 





where A(T) is the attenuation of the film at tem- 
perature T (normalized to the attenuation at T=0), 
and } is a small constant whose value depends 

on the theory. For the London theory b=1/6, 
while for the other theories b has not been eval- 
uated. Equation (1) is derived for the nonlocal 
theories on the assumption that the correlation 
length ~ is not dependent on temperature, and 

is valid if t<<,;. 

In the opposite limit where >>) 7 it can be 
shown that Eq. (1) is still valid if one substitutes 
Ae® for x7”. The quantity :,, is the penetration 
depth in the limit ~~. For the Pippard theory 
Ae “Az” so that the normalized attenuation is 
the same in both limits, but in the BCS theory 
Aw <rAz7/J(0, T), where J(0, 7) is a particular 
value of the kernel used in the BCS theory [see 
Eq. (5.45) of their paper]. Therefore the normal- 
ized attenuation is different in the two limits. 
The predictions of the BCS theory have not been 
evaluated for other values of é. 

Data were taken on four lead films varying in 
thickness by a factor 3, from 1.35°K to 4.2°K. 
The data were analyzed as follows. Equation (1) 
can be solved for \7(T)/x,(0), giving 


hp (TY/r, (0) =A-**[1 - ba(1-A)/2n,70)]- (2) 


Assuming that .7(7)/ (0) is the same for two 
films a and b, of different thicknesses d, and dp 
at the same temperature, we solved for B=bd,?/ 
2) 7(0) in terms of the experimentally known 
quantities A,, Ap, and dg/dp. We then calculated 


Q(T) =A~*[1 - B(1 - A)(@/d,)*], (3) 





53 








VoLUME 4, NUMBER 2 


PHYSICAL REVIEW LETTERS 


JANUARY 15, 1960 





for all four films at all temperatures. For this 
calculation one need only know B and the ratio of 
film thicknesses d/d,. [The quantities 6, »7(0), 
and d need not be known separately. If A;(0) is 
assumed to equal 400 A, then 5 is found to equal 
about 0.096. | 

If one assumes that B is independent of tem- 
perature, then 


Q(T) = A, (T)/a, (0) for E<<A, 


and 


Q(T) =A 7(T)/A7(0) for E>>a (4) 


L’ 
A plot of Q(T) calculated from the experimental 
data as discussed above is shown in Fig. 2. The 
fact that the data from all films fall on the same 
curve is an indication of internal consistency and 
attests to the validity of the method of analysis. 
The predictions of the Gorter-Casimir’’ two- 
fluid model and of the BCS theory for Az and 
d.”? are also shown in the figure. The BCS 
curves were calculated with the aid of tables 
compiled by Muhlschlegel.*® 

Curves of A,.2” and Ay were calculated for 
larger zero-degree energy gaps, assuming the 
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FIG, 2, Comparison of experimental curves of 
Q(T) vs y =(1-t*)~¥? with the predictions of Gorter- 
Casimir and BCS. The experimental curve should be 
compared with Ay in the limit £<<A;7 and with A”? in 
the limit £>>Ay. The reduced temperature ¢t =T/T,. 
was calculated for the experimental points using T. 
= 7. 22°K. 
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FIG. 3. Comparison of the experimental data with 
the BCS predictions of 47(7)/A;(0) for various energy 
gaps. 


temperature variation of the energy gap given by 


the BCS theory.’* The experimental data are con- 


pared with curves of A; in Fig. 3 and with curves 
of d,.2” in Fig. 4. Good agreement is obtained 
for an energy gap of 4.93kT in Fig. 3 and with 
a gap of 5.4k7,, in Fig. 4. The shape of the curve, 
however, indicates better agreement with the \ 
curves, which is not implausible in view of the 
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FIG. 4. Comparison of the experimental data with 
the BCS predictions of A,,”%(T)/ r,.2*(0) for various 
energy gaps. 
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fact that the correlation length é is limited by 

the thickness of the films. The energy gap found 
in this experiment is to be compared with the 
value of 4.127 ,, obtained from infrared absorption 
measurements.'® The apparent disagreement is 
not surprising, since our values were obtained 
using the BCS theory, which does not apply very 
well to lead. 

A detailed description of the experimental 
methods will be published in the April 1960 edition 
of the IBM Journal of Research and Development. 
We wish to thank Dr. A. G. Redfield and Mr. R. J. 
Blume for many helpful discussions and sugges- 
tions. 
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VARIATION OF SECONDARY ELECTRON EMISSION OF SINGLE CRYSTALS 
WITH ANGLE OF INCIDENCE* 


A. J. Dekker 
Institute of Technology, University of Minnesota, Minneapolis, Minnesota 
(Received November 23, 1959) 


Laponsky and Whetten’ recently reported re- 
sults of measurements of the secondary electron 
emission from single crystals of MgO and LiF 
as a function of the angle a between the incident 
electron beam and the normal to the surface. In 
previously reported work’ on polycrystalline 
materials the yield had been found to increase 
monotonically with increasing a; among other 
things, this is due to the larger production of in- 
ternal secondaries closer to the surface. The 
interesting feature of the results of Laponsky and 
Whetten lies in the occurrence of a number of 
maxima and minima superimposed on the “nor- 
mal” increase of the yield with a. The maxima 
occur at values of a for which the incident beam 
coincides with an important zone axis of the 
crystal. For MgO they find that the relative 
magnitude of the maximum corresponding to a 
=0 increases with increasing primary energy 
until it reaches a plateau of approximately 20% 
near 5 kev. 

Similar results have been obtained for titanium 
single crystals by Soshea® in our laboratory. 
However, he found that the relative magnitude of 
the maximum for a = 0 increases with increasing 


energy up to about 250 ev, and then decreases 
until it disappears near 12 kev. 

Since the angular position of the maxima is in- 
dependent of the wavelength of the incident elec- 
trons, Laponsky and Whetten conclude that the 
fine structure does not appear to be caused by 
electron diffraction. In the present Letter, we 
suggest that the effect does result from diffrac- 
tion of the incident beam. Since a comparison 
between the proposed theory and the experimenial 
results requires rather detailed numerical cal- 
culations, a more complete account will be given 
later together with the experimental results of 
Soshea. However, the essential features of the 
interpretation may be presented here. 

As the primary electrons penetrate the solid, 
they suffer inelastic collisions. Consequently, 
the depth over which coherent scattering can 
occur is limited to a surface layer with an effec- 
tive thickness of the order of the mean free path 
for inelastic scattering. The effective thickness 
thus constitutes only a fraction of the penetration 
depth of the incident electrons (which is ~300 A 
for a 1-kev electron in MgO). Presumably, the 
effective thickness increases with increasing en- 
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ergy in the energy range of interest (0.1 - 10 kev). 
Since for a small effective thickness of the crys- 
tal the von Laue condition corresponding to this 
direction is strongly relaxed, the incident beam 
will suffer diffraction of the same type as Kikuchi 
first observed in transmission experiments with 
very thin flakes of mica.* In such experiments 
the diffraction pattern consists of a number of 
spots on a film which are particularly intense 
and numerous whenever the direction of the in- 
cident beam coincides with an important zone 
axis of the crystal. Although the position of 
these spots depends on the energy of the incident 
beam, the effect as such is rather insensitive to 
this quantity. We believe that this kind of coher- 
ent scattering is responsible for the fine struc- 
ture mentioned above. 

In order to describe the diffraction of the in- 
cident beam, consider a lattice of identical 
atoms with primitive translation vectors a, b, 
¢. The positions of the atoms are then given by 
vectors r=n,a+n,b+m,c where n,, n,, n, are in- 
tegers. Let S, and 5 be unit vectors along the 
direction of the incident and scattered beams, 
respectively. The phase difference between the 
waves scattered by the atom at r and by the one 
at r=0 is then 


o= (20/r)G - 5,)*F = 2n(k - k,)*F =2n K-r. 


Here, \ is the wavelength of the electrons and 

K =2k,sin@; the scattering angle 6,=20. Let the 
surface of the crystal be parallel to the plane 
through a and b, and choose 1 cm? of the surface 
with edges N,a and N,b such that N, = N, = 107-10°. 
For a particular \ let the effective thickness 
along the ¢-axis be N,c. Furthermore, let K 
=h,a*+h,b*+h,c*, where a*, 6*, c* are the 
primitive reciprocal lattice vectors; h,, h,, h, 
are not necessarily integers. The intensity of 
the beam diffracted in the direction § is then 


E*(@)|>> exp(2niK-r) |? 


_ p29) [Sin?aN,h,| [sin?xN,h, | [sin?aNh 
Nt [ sin'sh, sin'sh, | = a (1) 





where E(@) is the atomic scattering factor. The 
functions in square brackets are each equal to 
N;? whenever h; is an integer (or zero); also, 
they vanish for the first time when hj =1/Nj. 
Thus, if N, is relatively small, the reciprocal 
lattice points are considerably extended in the 
direction of ¢*; along the extension, the factor 
in (1) containing N, decreases from N,’ to zero. 
The small subsidiary maxima occurring in these 
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functions between successive integer values of 
the h; may be neglected for our purpose. 

The effective thickness expressed by N, may be 
expected to be a function of \ and a. In the ab- 
sence of better information at present regarding 
the energy dependence of the mean free path for 
inelastic scattering it does not seem unreason- 
able to assume that N, is proportional to the 
primary energy E,, the proportionality factor 
being of the order of 10°? - 107° ev™*. Assuming 
one knows N,(d, @), the directions into which the 
primary beam is scattered for given \ and 5, are 
obtained as follows: Taking a particular reci- 
procal lattice point as origin, draw the vector 
-5,/r; the point so obtained is the center of the 
Ewald sphere of radius 1/.. The lines joining 
the center of the sphere with the intersections 
between the sphere and the extended lattice points 
correspond to the diffracted beams. This is in- 
dicated in Fig. 1 for a simple cubic lattice, neg- 
lecting the dependence of N, on a; we assumed 
N,=4X10°°E,, E,= 1000 ev, and a cube edge a=34 
Although in Fig. 1 only four layers of the recipro- 
cal lattice have been drawn, it is of course un- 
bounded. The fraction of electrons scattered into 
each of the diffraction directions can be calcu- 
lated from (1); for a simple cubic lattice one 
finds for this quantity approximately 


£*(8)(n2/a") cos(@, +a) ea (2) 


assuming that the incident beam lies in the (100) 





or=22.5° 


a=0° 


FIG. 1. Illustrating incident and diffracted beams 
as explained in text. 





1960 


of 


Ly be 


(2) 


00) 





VoLUME 4, NUMBER 2 


PHYSICAL REVIEW LETTERS 


January 15, 1960 





plane. 

In order to estimate the effect of the diffraction 
of the primary beam on the secondary emission 
we assumed that the diffracted beams would con- 
tribute to the yield in the same way as the “nor- 
mal” monotonic increase in yield with increasing 
angle of incidence predicts. The proposed theory 
explains qualitatively the experimentally observed 
features; calculations are in progress and have 
given satisfactory results so far. 

The author wishes to thank R. W. Soshea for 
stimulating discussions. 
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NUCLEAR RELAXATION VIA QUADRUPOLE COUPLING ms 


Robert L. Mieher 
Department of Physics, University of Illinois, Urbana, Illinois 
(Received November 30, 1959) 


Several papers'™ have reported calculations 
of the transition probabilities for quadrupolar 
nuclear spin-lattice relaxation and considerable 
data are becoming available.5-*° Nevertheless, 
three important problems concerning quadrupolar 
relaxation have not been solved. The first prob- 
lem is to express the relaxation time, 7,, in 
terms of the predicted transition probabilities in 
order to obtain a prediction that may be com- 
pared with experiment. The second problem in- 
volves the angular dependence of 7, as the 
orientation of the crystal is varied with respect 
to the magnetic field. Experimental 7,’s are 
isotropic, but the theoretical transition probabil- 
ities have a strong angular dependence. The 
third problem concerns the temperature depend - 
ence of the quadrupolar T,. Khutsishvili has 
proposed a high-temperature “four phonon” 
process as opposed to the Raman “two phonon” 
process considered by Van Kranendonk.’ This 
Letter is concerned with these three problems. 

Pound” first demonstrated and discussed the 
importance of the nuclear quadrupole moment 
in determining the spin-lattice relaxation time. 
The first detailed theory was developed by Van 
Kranendonk.’ Using the approximations that the 
lattice vibrations could be described by a Debye 
spectrum and that the lattice could be described 
by an array of point charges located at the lat- 
tice sites, he showed that the dominant quad - 
rupolar relaxation mechanism was a Raman 
process in which one quantum of lattice vibration 
is absorbed, another quantum emitted, and a 


nucleus makes a Am =:1 or :2 transition. He 
obtained explicit expressions for the Am =1 and 
2 transition probabilities, W, and W,, for the 
NaCl lattice. 

The spin systems we consider are those in 
crystals of cubic symmetry with no static quad- 
rupole interactions. The nuclear energy levels 
are equally spaced and are given by the Zeeman 
interaction, E(m)=-yhHm. Abragam and Proc- 
tor*® have shown that the concept of spin tem- 
perature is valid for such systems when the 
interactions between nuclei are stronger than 
the interactions between the lattice and the nuclei 
(i.e., T,“«T7,). Using the spin temperature con- 
cept, an explicit relation between relaxation 
times and transition probabilities has recently 
been given by Hebel and Slichter.* This relation 
is 


Do. (E_-E 


2 
a m+p) Vim, m+n) 
772 SE? » & 

mm 





where the energy values, E,,, are chosen such 
that the trace of the energy matrix is zero and 
W(m, m+) is the probability of a nucleus making 
a transition from state m to m+u. 

Starting with the quadrupole Hamiltonian of 
Cohen and Reif,'® we have carried out calcula- 
tions for the transition probabilities for a point- 
charge zincblende lattice. The predicted value 
of W(m, m+) is 


W(m, m+) =C\ a PTME(T)E (=), (2) 
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where C =V3/8nd*y°r*, r = interatomic distance, 
d = density, v = velocity of sound, T* = T/@ 

(@ is the Debye temperature), Q,,,, is the matrix 
element for the quadrupole operator between 
states m and m+u, E(T* ) is a numerical func- 
tion given by Van Kranendonk, and the E pl@) 

are 


™ (rey 2 
Ej (=) =“4e- [723.4 - 3120"), 


(vey 2 
E (@) =" =~“ [645.4 + 787], (3) 


where a? =a,7a,7+@,7a," + 0,70,” and Q,, Q2, 
are the direction cosines between the magnetic 
field and the [100], [010], [001] cyrstalline direc- 
tions. The notation is the same as that used by 
Van Kranendonk. 

When the expressions for the transition prob- 
abilities are substituted into Eq. (1), the result 
is 


F780 FOP 1) OCT ETE (=) +4851), (4) 
where the term / (J) = (2/+ 3)//*(21 -1) contains 
the entire dependence of 7, on the nuclear spin. 
This spin function is generally valid for the 
NaCl, CsCl, and ZnS lattices because the lattice 
details are contained in the E p() terms. While 
it has been recognized that 7, would depend on 
the nuclear spin,® the explicit dependence has 
not previously been derived. 

When the values of £ ,,(~) are inserted in Eq. (4), 
it becomes 


1 1 (2%+3) (e7)@)? 


T, 40 7(27-1) r° 

T, is isotropic because the angular dependences 
of the transition probabilities exactly cancel. 
This isotropy holds also if one uses Van Kranen- 
donk’s* transition probabilities for the NaCl 
lattice or Yosida and Moriya’s® transition prob- 
abilities for the NaCl and CsCl lattices. [Be- 
cause of a different choice of quadrupole Hamil- 
tonian their coefficients of a? in E,(~) are 1/16 
those of a” in E,(~). This is balanced by a factor 
of two difference in the Qum’s so that the pre- 
dicted relaxation is still isotropic.] Our T, 
measurements on In" in a single crystal of 
InSb at 77°K give an isotropic 7, of 0.80+0.05 
second. Other investigations®»® have reported 
isotropic 7,’s in the alkali halides. 

The spin function, f(/), has been experimentally 
verified with measurements on AlSb and RbCl. 


CT**E(T*)(3305). (5) 
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These results are 





Ce 3, ye 0.425, 
ris) : (6) 7 = 0.69 (theory), 
raise 7.75 + 0.10 (exp.); 
[eer -a.2e, £82) 0.24, 
T,(Rb*") 


. 4 =1.027 (theory), 


T,(Rb®) f 


qT, (Rb®*’Cl) 
T,(Rb*C1) 
Figure 1 shows the results of temperature 

dependence measurements of I'*” in KI between 
77°K and 800°K. The data lie close to the 
T,«T™ line, the predicted high-temperature 
dependence for the Raman “two phonon” process. 
The 7, «7™ line is arbitrarily drawn through 


=1.23 + 0.40 (exp.). 
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FIG. 1. Relaxation times of I?’ in KI. The “two 


phonon” process predicts T,« 7~ and the “four 
phonon” process predicts 7, « T* 
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FIG. 2. Relaxation times of In‘ in InP. Because 
T<6, the data deviate from the straight line, T,;< 7~’; 
the dashed line is calculated using Van Kranendonk’s 
function, E(T*). 


300°K. It is apparent that the 7, dependence of 
the “four phonon” process proposed by Khutsish- 
vili’ does not fit the data. 

The temperature dependence for the Raman 
process deviates from T~ for 7 < 6 and is given 
by Van Kranendonk’s function E(T"). Figure 2 
shows the results of 7, measurements on In*** 
in InP. The straight line is T, «T~* and the 
dashed line is the predicted temperature de- 
pendence using E(T*). The line is calculated 
using the room temperature 7, and a Debye 
temperature of 400°K which gives the best fit to 


the data. The scatter in the data gives a + 25°K 
uncertainty for the Debye temperature (which is 
not available from other sources). 

The author wishes to thank Professor C. P. 
Slichter for numerous suggestions. Some of the 
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PRESSURE DEPENDENCE OF THE CURRENT-VOLTAGE CHARACTERISTICS OF ESAKI DIODES 


«7 


S. L. Miller and Marshall I. Nathan 
International Business Machines Research Laboratory, Poughkeepsie, New York 


and 


Arthur C. Smith 
Department of Electrical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 
(Received December 8, 1959) 


Esaki’ has discovered a new tunneling phenom- 
enon in narrow p-n junctions which yields a neg- 
ative resistance characteristic in the forward 
direction. Measurements are reported here on 
the variation of this phenomenon with hydrostatic 
pressure in germanium. These data make it pos- 
sible to obtain the dependence of tunneling current 
on the energy gap E g Measurements of the “or- 
dinary diode current” in the forward direction 
indicate that the pressure variation of E, is 
approximately the same for degenerate material 
and pure germanium. 

Pressure experiments at room temperature 
were carried out in an apparatus of the type de- 
scribed by Bridgman’ using pentane as a pres- 
sure transmitting fluid. 

The variation of the peak current with pres- 
sure divides all of the diodes studied into three 
distinct groups. Representative curves for these 
groups are shown in Fig. 1 in which InJp is 
plotted against pressure. It is found that almost 
the same curves result from a plot of (dJ/dV) y_ 0 
vs pressure. Those diodes whose J vs V curve 
at 4°K exhibits structure caused by phonon par- 
ticipation in tunneling® give curves like (a). These 
are very close to straight lines. Curves similar 
to (b) and (c) are obtained for units in which there 
is no evidence for phonon-assisted tunneling at 
4°K. Curves of type (c) are typical of units made 
by alloying n-type material (SnAs) on a p-type 
substrate while the units giving curves (a) and 
(b) were made by alloying p material (either SnGa 
or InGa) on an n-type substrate. The total rela- 
tive change of peak current with pressure in- 
creases with decreasing current density from 
diode to diode. It varies from a factor of 5 to 
about 65 in 30000 kg/cm?. The break in curves 
of type (b) occurs at 8000 kg/cm? regardless of 
the magnitude of J, (over at least a factor of 
seven), the kind of impurity in the substrate, or 
the nature of the alloying dot. The pressure 
effects are independent of the crystallographic 
orientation of the substrate. The insensitivity 
to orientation and the results of a resistivity 
measurement on a sample of the same geometry 
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lead to the belief that the pressure is essentially 
hydrostatic. 

The valley current, J, (Fig. 1), is larger than 
is predicted on the basis of a simple tunneling 
picture and ordinary diode theory.* J, also 
shows a decrease with pressure, although the 
over-all change is smaller than the change in 
Jy for the same unit. There is frequently a plot 
like curve (b) with a break at 8000 kg/cm?. This 
is, however, not correlated with the detection 
of phonon-assisted tunneling current. It may be 
an indication that the valley current is related 
in physical mechanism to the nonphonon tunneling 
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FIG. 1. Peak current density divided by its value at 
atmospheric pressure vs pressure. The circles, 
triangles, and squares designate points taken on the 
return to atmospheric pressure. Inset shows typical 
J -V characteristic of these diodes. 
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current. 

In the region of the characteristic where the 
“ordinary diode forward current” predominates, 
it is found that J« exp(eV /kT) at room temper - 
ature. If the usual diode equation is assumed, 
the change in V at constant current will equal 
the change in E g with pressure provided that 
variations in effective mass, lifetime, and mo- 
bility are neglected.® The variation in E, with 
pressure up to 30000 kg/cm? obtained by this 
method is the same to within 15% for these 
heavily doped materials as that observed in 
intrinsic Ge by measurement of bulk conductiv- 
ity. These data are plotted in Fig. 2 for two 
diodes. Since the effective energy gap for this 
phenomenon is determined by the energies and 
effective masses of the (111) and (100) conduction 
band minima, this is an indication that the band 
structure is not seriously perturbed. 

Some features of the data in Fig. 1 are con- 
sistent with the tunneling theory given by Keldysh’ 
and Price and Radcliffe.* Their results for 
phonon-assisted tunneling can be put in the form: 


J, =Kexpl-ak W,)- (1) 


W,= junction width at atmospheric pressure and 
zero applied bias. a is a parameter which de- 
pends on the effective mass and the dielectric 
constant. In what follows, the pressure varia- 
tions of a and K are neglected despite the fact 
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FIG. 2. The change in voltage with pressure at con- 
stant current compared with the known change in energy 
gap with pressure for intrinsic germanium. 


that these may occur. From Eq. (1) it can be 
shown that 


InJ A, Ink +E,(dind ,/dP)/ (dE,/aP), (2) : 


where the derivatives are with respect to pres- 
sure. 

Fig. 3 shows a plot of Ind, vs (dinJp/dP) for 
all of the units measured in the pressure experi- 
ments except those made by alloying ” on p ma- 
terial. A line drawn through the points for the 
units showing phonon-assisted tunneling has a 
slope of 12.4x10* kg/cm?. This number agrees 
with the value obtained by taking 


E,/(dE,/aP) = 0.65 ev/(0.05 x 107* evcm?/kg) 
=13x10* kg/cm’. 


A determination of W, for séveral selected units 
by measurement of capacitance leads to a value 
for the parameter a in Eq. (1) of 2x10" (ev- cm) 
for phonon-assisted current. This corresponds 
to a reduced effective mass of 1/20 of the free 
electron mass, if Keldysh’s’ form for a is used, 
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FIG. 3. Peak current density at atmospheric pres- 
sure vs (dind p/4P). 
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and a field of E,,/W, is assumed. 
It can easily shown that if the phenomenol- 
ogical law, 


_~ nr 
J, =Kexp-f5 | ), (3) 


governs tunneling, then the slope of the plot in 
Fig. 3 would be divided by v. It is apparent that 
the data on phonon-assisted tunneling obviate 
values of nm which depart greatly from unity. On 
the other hand, the plot for the low-pressure 
slope of units (b) seems to fit with Eq. (3) for 
n=2. This plot is also consistent with Eq. (2) if 
Eg is interpreted to be the gap between the (000) 
minima of the valence and conduction bands.° 
This suggests that the process involved may be 
a second order one in which there is virtual 
tunneling between the zone center states.'° 

Measurements on two units of type (b) at liquid 
N, temperature to a maximum pressure of 5000 
kg/cm? indicate that the pressure dependence of 
InJp is about 25% larger at the low temperature 
than at room temperature. This is not easily 
understood on the basis of a phenomenological 
law like Eq. (2), since Jp is almost independent 
of temperature. However, these low-temperature 
pressure measurements are complicated by the 
possibility of nonhydrostatic effects. More care- 
ful measurements are necessary. 

There are two general features of the data in 
Fig. 1 which are not understood at the present 
time. The first is the abrupt change in slope of 
units (b) at 8000 kg/cm?. Nothing drastic is known 
to happen to the band structure at this pressure. 


The second is the difference in shape of curves 
(c) and (b). 

We wish to express our appreciation for the 
hospitality extended to us by Dean Harvey Brooks 
and the members of the high-pressure group at 
the Gordon McKay Laboratory, Harvard Univer- 
sity. We thank our colleagues at the IBM Re- 
search Laboratories for stimulating discussions 
and for assistance in fabricating diodes. 
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SPIN RESONANCE OF CONDUCTION ELECTRONS IN InSb 


G. Bemski 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received December 2, 1959) 


Among the interesting properties of InSb is 
the very high and negative magnetic moment of 
its conduction electrons.' 

We have observed directly spin resonance of 
conduction electrons in n-type InSb in the con- 
centration range between 2 x10“ electrons/cm* 
and 3x 10'5 electrons/cm*. The |givalues varied 
between 48.8 and 50.7 depending on the electron 
concentration. The experiments were performed 
at temperatures of 1.2°K and 4°K. By observing 
the resonance at two different frequencies (9000 
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Mc/sec and 24000 Mc/sec) we verified that it 
was due to spin rather than plasma resonance. 
The magnetic field at which the resonance was 
observed was directly proportional to the micro- 
wave frequency while the plasma resonance gives 
rise to an inverse dependence of the magnetic 
field on the electronic frequency.? Cyclotron 
resonance can also be ruled out on the basis that 
the observed resonance occurs at approximately 
2.5 times the field value of the cyclotron reso- 
nance. The widths of the observed lines (see 
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Table I. Experimental and calculated results. 








Electron Fermi 

concentration level 
(em~) (ev) Il exp (m*/M) ox, “(m"/m»)theor® AH 
2x 10" 9.6 10~ ev 50.7 0.0131 0.0130 29 
7x 10" 2.2107 ev 50.3 0.0132 3¢ 
1.510% 3.6 107 ev 49.6 0.0134 5@ 
2.510% 4.91078 ev 48.8 0.0136 0.0136 8” 
3x 105 5.5 107 ev 48.8 0.0136 23¢ 





8gee reference 7. 


Table I) are also more than an order of magnitude 
less than in the case of cyclotron resonance.* 

Figure 1 shows the barely resolved spin reso- 
nance line superposed on a plasma line. Figure 
2 shows the spin resonance line greatly amplified. 
In this case the sample, 10-? cm* in volume, was 
positioned directly at the end of the waveguide. 
Runs were also performed with much smaller 
samples in rectangular cavities. The lines are 
isotropic with respect to the field orientation. 

The InSb band structure has been treated by 
Kane who has derived expressions for curvature 
of the nonparabolic conduction band.* Cyclotron 
resonance has shown that the conduction electrons 
possess a low effective mass, m*=0.013m,, near 
the bottom of the conduction band.* A variety of 
other experiments have given values of m* as a 
function of the energy of electrons above the 
bottom of the band.® 


eo 


; ; ; ; ‘ ; i 
4 oe oe | $+ woe ye +. 
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FIG. 1. Plasma resonance line in ”-type InSb at 
1.2°K. Encircled is the superposed electron spin 
resonance line. 





Roth et al.’ have calculated the expression for 
the magnetic moment pz: 


__eh (m™ _\._4 
mes \m* yEasT\t 


where A is 





the spin-orbit interaction=0.9 ev as 


(a) 








FIG. 2. (a) Electron spin resonance absorption in 
n-type InSb with 2 x 10" electrons/cm*. (b) Similar 
resonance absorption in a sample with 3 x 10" electrons 


per cm, 
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calculated by Kane, and E g is the energy gap of 
0.25 ev at helium temperature.® Zwerdling, Lax. 
and Roth’»® have observed magneto-optical transi- 
tions to the lowest Landau Level (/=0) in the 
conduction band which gave the value of p =26p p.* 
In the present experiments the shift of the g 
value with concentration is also believed to be 
associated with the increasing mass of the elec- 
trons. Table I shows the close agreement be- 
tween the observed g’s and those calculated by 
using Eq. (1). The values of m* were calculated 
from the known Fermi energy in the samples.’ 
All samples were degenerate at helium temper- 
atures with the purest sample having a degeneracy 
temperature of 11°K. Since the only electrons 
which participate in the resonance absorption 
belong to the highest occupied Landau level, the 
amplitude of the signal does not increase with the 
concentration of electrons. We were unable 
to observe resonance in a sample with 10** 
electrons/cm*. This is probably due to the pro- 


gressive broadening of the lines with concen- 
tration (Table I). 

I would like to thank C. Lee and R. Logan for 
the InSb samples, G. Feher for the use of the 
X-band spectrometer and B. Szymanski and E. 
Gere for help with the experiments. 
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MILLIMETER CYCLOTRON RESONANCE IN SILICON* 


C. J. Rauch, J. J. Stickler, H. J. Zeiger, and G. S. Heller 
Lincoln Laboratory, Massachusetts Institute of Technology, Lexington, Massachusetts 
(Received December 7, 1959) 


A cavity spectrometer operating at 2.1 mm 
(136 kMc/sec) has been constructed and cyclo- 
tron resonance absorption experiments have 
been performed on high-purity n-type silicon 
from 1.2°K to about 50°K. 

About 50 microwatts of 2.1-mm microwave 
energy are generated as the second harmonic 
from a crossed guide crystal harmonic genera- 
tor which is fed with 4-mm fundamental power 
from a Philips DX 151 klystron. The harmonic 
output is transmitted through an oversized wave- 
guide to a circular 7M,,, cavity located at the 
bottom of a helium Dewar. The sample, which 
is about 0.020 in. diameter and about 0.040 in. 
long, is coupled to the high electric field region 
of the cavity through a small axial hole beyond 
cutoff in the end wall of the cavity. The sample 
is fastened to a quartz rod which extends to the 
top of the Dewar and serves to position the sam- 
ple in the cavity and as a light pipe for carrier 
excitation. Temperature is monitored with a 
carbon resistance thermometer fastened to the 
top of the cavity. The cavity, waveguide, and 
quartz rod are encased in a shield to prevent 
direct contact with liquid helium, thus eliminat- 


64 


ing noise due to bubbling of liquid helium in the 
system. A low pressure of He gas inside the 
shield provides heat transfer from the cavity 
and sample to the helium bath. Stabilization of 
the klystron frequency is accomplished by apply- 
ing a small 1-kc/sec signal to the reflector of 
the klystron, and feeding back the synchronously 
detected signal to the reflector. Resonance lines 
are observed as a function of the static magnetic 
field in the standard way by modulating the car- 
rier density with a light chopper at a low audio 
frequency and synchronously detecting the signal 
as the magnetic field is swept. 

Very pure n-type silicon (3000 ohm cm at 
room temperature) has been used to obtain the 
effective-mass parameters of the constant ener- 
gy ellipsoids associated with the conduction 
band minima. Orientation data were taken to 
obtain experimental values for the sum 1/m,? 
+1/m,? +1/m,’ (see Fig. 1) which is constant for 
all orientations and is given by 
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ported values; however, our m7; lies outside the 














baa | ; limits of error of the values previously reported.’ 
i | “ga Ware wa. ae At present no explanation of this discrepancy 
a | a pe, can be given. 

7 Measurements of the effective mass and line 


width were also made as a function of tempera- 
ture from 1.2°K to 50°K with the crystal oriented 
so that the static magnetic field was approxi- 
mately along the [100] direction. No indication 
of a change in mass was observed. wt asa 
function of temperature measured at low light 
intensities is shown in Fig. 2. (wr =2H,/AH, 
where H, is the field at resonance and AH is the 
line width.) Values of w7 approaching 200 were 
observed, corresponding to a “cyclotron mobility” 
of the order of 1.5x10° cm? /volt-sec. 

Also shown in Fig. 2 is a curve of m*wu/e, 
where m* is given by 1/m”™ = 3(1/m)+2/m;}) and 
pu =4.0x10°T -**® cm? /volt-sec is the high-tem- 

a —— — -— = a $0 perature conductivity mobility for conduction 

RELATIVE CRYSTAL ROTATION ANGLE (degrees) electrons given by Morin and Maita.? Our ex- 
perimental results could be extrapolated to the 
FIG. 1. Experimental values of 1/m,*+ 1/m,?+1/m,?, straight line m*wu/e. However, we were un- 
2(1/m,’+1/m,"), and 1/m;” (in normalized mass units) able to obtain data at temperatures higher than 
as a function of sample orientation at 4. 2°K. 50°K due to the masking of the light-excited car - 
riers by intrinsic carriers. 

it is hoped that the improved resolution (about 

6 times that of a K-band system) and high sen- 








Here m,, m2, and m, are the effective masses 
measured for the individual ellipsoids in the 
conduction band of Si; m; and m, are the trans- 






































verse and longitudinal components of the ellip- ec! RR | TTT wee 
soids, and m* [111 ] is the effective mass in the | \| e 
200 ttt t | 
[111] ‘nent The maximum of the sum Te | . on | 
3[1/m,? +1/m,”] as a function of orientation me mt Mer | 
yields the effective mass m *[110] given by =, 1 | 
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In addition, an extra relation is found by plotting “ | Al YI] | 
the experimental values of 3(1/m,?+1/m,”) and . — ; a mn ts eS 
1/m,’ as a function of orientation. Where the two ° -—- rt | 
curves intersect, the following holds: td — : y = 
| 
1 1 1/1 || |_| \ | 
bac —_ 2) Cram, m } (3) 2] T t | T *% | 
3 Jintersection [111] t til | | Bai ge ae Pilia ee 
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Figure 1 shows experimental curves of 1/m,’ we 
+1 /m? +1 /m2, (1 /m? +1 /m.?), and 1 /m;? as FIG. 2. wrt as a function of temperature for cyclo- 
a function of sample orientation at 4.2°K. Using tron resonance absorption in high-purity n-type silicon 


at 135.95 kMc/sec. The observations were made on 
: ; the low-field electron resonance with the magnetic 
with Eqs. (1), (2), and (3), the self-consistent field approximately along the [100] axis. (The dif- 
values mz =(0.192 + 0.001)m, and m;=(0.90+0.02)m, ferent symbols refer to different runs.) Also shown 
are obtained, where m, is the free electron mass. is a curve of m*wy/e, where p is the high-tempera- 


Our value of m; agrees with previously re- ture conductivity mobility for electrons. 


the appropriate values from the curves in Fig. 1 
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sitivity of the 2-mm spectrometer will permit 
the observation of cyclotron resonance in other 
materials. 

We wish to thank Dr. Benjamin Lax for his 
continued interest in this work. 
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THRESHOLD ENERGY FOR LATTICE DISPLACEMENT IN a-Al,0O, 


George W. Arnold and W. Dale Compton 
United States Naval Research Laboratory, Washington, D. C. 
(Received December 7, 1959) 


Seitz' has estimated that the threshold energy, 
Ey, for displacement of an ion or atom ina 
tightly bound solid should be about 25 ev. Ex- 
periments have been carried out on metals*~® 
and semiconductors®’ which yield displacement 
energies of 10 to 30 ev. This Letter reports 
the first direct determination of £ for a non- 
conducting material, i.e., a-Al,O,. 

Levy and Dienes® have shown that reactor 
irradiation of a-Al,O, causes the growth of an 
optical absorption band at about 204 my which 
cannot be generated by x-rays or » irradiation. 
This band is presumed to arise from the trap- 
ping of charges at defects which occur when 
atoms are displaced. 

In the present experiment, 0.0635-cm thick 
samples of a-Al,O, (Linde) were irradiated at 
77°K with fast electrons from a Van de Graaff 
accelerator and the resulting optical absorption 
was measured at 77°K without intervening warm- 
up with a Cary Model 14M spectrophotometer. 
In Fig. 1 are shown the absorption features in- 
duced by irradiation with 2-Mev electrons at 
77°K and by a fast-neutron irradiation® at pile 
temperature (Oak Ridge X-10 reactor). As 
previously reported by Levy and Dienes,® ab- 
sorption bands were found at 205, 230, and 255 
my in neutron-irradiated material. The only 
band observed after electron irradiation was at 
205 mu. The absorption on the long-wavelength 
side of the 205-my band in the electron-irrad- 
iated sample bleached when the sample was 
warmed to room temperature, and a slight in- 
crease in absorption took place at 205 mu. No 
bleaching of the 205-my absorption occurred 
when the samples were irradiated and held at 
77K. A very slow bleaching of this band took 
place when held at room temperature after the 
low-temperature irradiation. Irradiation at 
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77°K produces many more centers than does a 
comparable irradiation at room temperature, 
the ratio of the yields at these two temperatures 
being at least ten. 

The growth of the 205-myz band with increas- 
ing electron flux was qualitatively similar to 
that found by Levy and Dienes® for neutron ir- 
radiation, i.e., a rapid initial increase followed 
by a linear growth. No dependence on rate or 
crystallographic orientation with respect to the 
electron beam was found. 

The formation of the 205-myz centers as a 
function of electron energy was determined at 
a constant incident flux (0.945 x10" e/cm?). 
This flux gives absorption values on the rapidly 
rising portion of the growth curve. The energy 
absorbed in the sample, £’, is obtained from 
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FIG. 1. Optical density vs wavelength in my for 
a@-Al,O;. Background prior to irradiation has been 
subtracted. —— 10!" mvt (fast), pile temperature. 

--- 2-Mev electrons, 9.45 x 10'* e/cem?, 77°K. Meas- 
urements made at 77°K. 
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the stopping power of a material’® of effective 
atomic number Z =11 and from multiple scatter- 
ing theory."* Until the energy of the slowed- 
down electron falls below the threshold energy, 
displacements can occur. 

The effective energy, E eff’ is the mean energy 
of the electron in decreasing from its value at 
the surface to the threshold energy. If the elec- 
tron emerges from the back surface with an 
energy greater than the threshold energy, E ors 
is the mean energy it possesses while traversing 
the sample. A quantity a is defined as the meas- 
ured optical density (corrected for background 
and overlap of other bands) divided by the cor- 
rected path length. This quantity should be in- 
terpreted as the number of centers formed per 
cm? per unit path length of the electron while 
it produces defects. This path length is de- 
termined as follows: (1) If the energy of the 
electron exceeds the threshold, the path length 
is determined from the theory of multiple scat- 
tering.’ (2) If the energy falls below the thres- 
hold, the path length is computed from the stop- 
ping power of the material.*° In the latter case, 
small errors in determination of the threshold 
energy are magnified in computing this distance. 
No data at energies lower than 0.75 Mev are in- 
cluded because of the increasing error in this 
computation as the incident electron energy 
approaches the threshold value. These quanti- 
ties are tabulated in Table I. 

The ratio of the quantity a at the various ef- 
fective energies to that at 1.77 Mev is plotted 
in Fig. 2 and compared with the calculated total 
cross sections, all normalized to unity at 1.77 
Mev, for aluminum and oxygen for various 
values of Eg.” The data are best fitted by the 
theoretical curve for E, = 50(+5) ev in the case 
of aluminum and by Ez =90(+5) ev in the case of 
oxygen, giving a threshold energy of 0.43 Mev. 


The data do not allow an unequivocal selection 
of the defect giving rise to the 205-my absorp- 
tion as being generated by either aluminum or 
oxygen displacements. Aluminum displacements 
might be more likely because they have greater | 
freedom of movement that do the close packed 
oxygen atoms. 

Table II gives the oscillator strength of the 
center that is obtained, using either aluminum 
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FIG. 2. Optical density per unit path length vs 
effective electron energy, compared with calculated 
total cross sections for oxygen and aluminum, all 
normalized to unity at 1.77 Mev. #, optical density 
per unit path length induced by a constant flux of 
9.45 x 10'6 e/cm? for various effective electron ener- 
gies. —— Calculated total cross sections for oxygen 
and aluminum for various values of Ej. 


Table I. The quantity a, normalized to unity at E,¢¢=1.77 Mev, for various incident electron energies. 














E 
Incident E’ Est 
electron Energy absorbed Effective energy 
energy in the sample of electron® 
(Mev) (Mev) (Mev) a aE efs)/Q(E ogg = 1. 77) 
2.00 0.46 1.77 3.14 1 
1.50 0.49 1.25 2.28 0.73 
1.00 0.59 0.72 1.08 0. 34 
0.75 0.73 0.59 0.81 0.26 





*Effective energy Eorp=E -E'/2 or Eefg=(E + 0.43)/2, whichever is larger. 
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Table Il. The oscillator strength of the defect giving 
rise to the 205-my absorption band, assuming that the 
defect arises from either an aluminum or oxygen dis- 
placement. 








Ed Oo” Nd 
Element (ev) (barns) (cm~*x10%) fF 
Aluminum 50 12.6 5.59 0.51 
Oxygen 90 4.20 2.80 1.02 





or oxygen as the defect in question. The quantity 
fNq is obtained from the optical measurement by 
using Smakula’s equation, where f is the oscilla- 
tor strength and Nj is the number of defects /cm’. 
Nq is obtained from Nqg=N,.0q?, where N, is the 
number of atoms/cm’, oq is the total cross sec- 
tion for displacement,” and ¢ is the incident 
electron flux ine/cm?. The values of og and Nq 
were calculated for an effective energy of 2.0 
Mev. The contribution of secondary displace- 
ments to Ny has not been included and is esti- 
mated to be less than 10% of the total. 

Dienes and Vineyard** have computed the 
oscillator strength of the 205-my absorption 
bands in neutron-irradiated a-Al,O, by a simi- 
lar calculation. This is based on the measure- 
ments of Levy and Dienes.*® Using a value of 
E4=25 ev they get f =0.011. A comparison of 
this oscillator strength with that determined in 
other materials led them to conclude that this 
value was unreasonably small. The irradiations 
were made at pile temperature. According to 
the results of the present experiment, only about 
one tenth of the centers generated at this tem- 
perature are retained. A re-evaluation of their 
experiment using a tenfold increase in the num- 
ber of centers, E,=50 ev, and considering only 
aluminum displacements, gives a value of f =0.34. 
This is in good agreement with the present ex- 
periment, especially since the fast-neutron flux 
value can be in error by 50%.'* 

Two other experiments have been made on 





- 


a-Al,O, in which comparisons can be made with 
collision theory. They are: (a) the experiment 
of Antal and Goland™ on the attenuation of a 
thermal neutron beam in passing through a fast- 
neutron-irradiated sample, and (b) that of 
Martin,*® who measured the external changes in 
dimensions of irradiated material. In both of 
these experiments fewer defects are obtained 
experimentally than are predicted by theory; 

in (a) by at least a factor of 40 and in (b) bya 
factor of 50. Modification of their results by 
the temperature factor and using Ej =50 ev 
brings experiment and theory into reasonable 
agreement. 
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PRECISION DETERMINATION OF THE HYPERFINE STRUCTURE OF THE GROUND STATE OF 
ATOMIC HYDROGEN, DEUTERIUM, AND TRITIUM* 


L. Wilmer Andersont and Francis M. Pipkin? 
Lyman Laboratory, Harvard University, Cambridge, Massachusetts 


and 


James C. Baird, Jr. 
Mallinckrodt Laboratory, Harvard University, Cambridge, Massachusetts 
(Received November 12, 1959) 


The hyperfine splittings of the hydrogen iso- 
topes have been measured by an optical polar- 
ization method.'~* In this experiment a small 
quantity of one of the hydrogen isotopes was 
placed in a 500-cm® spherical flask together 
with a buffer gas and rubidium metal. Rubidium 
resonance radiation produced in an electrodeless 
radio-frequency light source was circularly 
polarized, passed through an interference filter 
to select the D, line, and then used to polarize 
the rubidium. A pulsed 30-Mc/sec oscillator 
was used to dissociate some of the hydrogen 
molecules into atomic hydrogen. This atomic 
hydrogen became oriented by spin exchange col- 
lisions with the polarized rubidium. The polar- 
ization of the rubidium was monitored by observ- 
ing the transmitted resonance radiation with a 
935 phototube. When a radio-frequency field 
was applied so as to destroy part of the hydrogen 
polarization, there was a decrease in the amount 
of light incident upon the photocell. 

The following system was used to generate and 
measure precisely the frequencies for observing 
various transitions. The 1-Mc/sec output of an 
Atomichron, operated and standardized by J. A. 
Pierce of the Cruft Laboratory, was multiplied 
in a Gertsch AM-1 and used to phase lock the 
Gertsch FM-4 at a high frequency. The internal 
low-frequency oscillator of the AM-1 was not 
used. A T-116/APT-5A radar jamming oscil- 
lator was used to drive the hyperfine transitions. 
An auxiliary output of this oscillator was beat 
with the output of the Gertsch FM-4 obtaining a 
beat frequency between 3 and 4 Mc/sec. This 
frequency was amplified and compared in a phase- 
detecting circuit with a variable frequency Gen- 
eral Radio 616D oscillator, which was very 
Stable. This beat frequency was measured with 
a Northeastern Engineering 14-20 electronic 
counter, whose 100-kc/sec time base was sup- 
plied by the Atomichron. The output of the phase 
detection circuit was used to phase lock the APT - 
5. A relay was used to square-wave modulate the 


output of the APT-5 which induced the radio- 
frequency transitions. The output of the photo- 
cell was amplified and fed to a lock-in detector 
with an output time constant of 0.1 sec. 

The field-independent 0-0 transition in hydro- 
gen and tritium could not be detected in this 
fashion. This is probably due to the manner in 
which the spin exchange collisions populate the 
magnetic sublevels. For these isotopes it was 
necessary to measure the two field-dependent 
hyperfine transitions (F =0, m =0-F=1, m=+1) 
and to average the results to obtain a field- 
independent quantity. A typical line profile for 
hydrogen is shown in Fig. 1. Each of the field- 
dependent hyperfine lines is asymmetrical; this 
asymmetry, however, occurs on the high-fre- 
quency side for one line and on the low-frequency 
side for the other line. It is assumed in this 
experiment that when the frequencies of the two 
lines are averaged, the effect of the asymmetry 
cancels out. All of the measurements on these 
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FIG. 1. Typical line profile for the (F=0, m=0 
— F=1, m=1) hyperfine transition of atomic hydro- 
gen in a molecular hydrogen buffer gas. The lock- 
in-detector time constant was 0.1 sec. 
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two isotopes were made between the hours of 1 
and 5 a.m. when the magnetic activity (i.e., 
trolleys) in the neighborhood was a minimum. 
In deuterium, however, there are two transi- 
tions (F =3/2, m=1/2~-F=1/2, m=-1/2) and 
(F =3/2, m=-1/2~F=1/2, m=1/2) which are 
almost field independent (Fig. 2). These two tran- 
sitions differ by twice the product of the nuclear 
moment and the field. These two lines could 
easily be resolved in the field of 0.15 gauss which 
was used in the experiment and had a linewidth of 
about 175 cycles/sec. This residual width was 
produced by exchange collisions between the 
deuterium atoms. To determine the hyperfine 
splitting of deuterium, both these lines were 
measured and the Zeeman transitions (AF = 0) 
were used to compute the second order correc- 
tion. 

Since the value of the hyperfine splitting ob- 
served depends upon the pressure and nature 
of the buffer gas used, several bulbs were made 
at various pressures for each gas. Each of these 
bulbs was then measured on a number of differ- 
ent nights. A least-squares fit was used to deter- 
mine the best straight line through these points. 
In this manner the pressure shift for each gas 
could be determined and the measured hyperfine 
splitting could be extrapolated to zero pressure. 
Table I shows the pressure shift obtained for 
the various buffer gases. Table II shows the 
values of the hyperfine splitting obtained by 
extrapolating to zero pressure. In the case of 
deuterium, where the observed transitions are 
field independent, the errors appear to be purely 
statistical. The value obtained for the zero-field 


Table I. The pressure shift data. 


— 327,384696 
— 327.384521 


FIG. 2. Typical line profile for the (F =1/2, m=-1/2 
—F =3/2, m=1/2) and (F =1/2, m=1/2—F =3/2, m = -1/2) 
hyperfine transitions of atomic deuterium in a neon : 
buffer gas. The separation of the two peaks is equal 
to twice the nuclear resonance frequency for the 
deuteron in a field of 0.15 gauss. The markers were 
made with only one direction of frequency sweep and 
serve only to mark the separation of the two peaks and 
the linewidth. The lock-in-detector time constant 
was 0.1 sec. 


hyperfine splitting is 
Av (H?) = 327.384349 + 0.000005 Mc/sec. 


In the case of hydrogen and tritium, where field- 
dependent lines are measured, the errors do not 
appear to be statistical. The amount by which 


The pressure in each bulb was measured at room temperature with an oil 


manometer when the bulb was prepared; the measurement of the hyperfine splittings were made with the bulb 
at 50°C. In computing the pressure shifts at 50°C the perfect-gas law was used to correct the pressure. The 
error in the pressure shift is estimated to be about 20%. 

















Hydrogen Pressure shift Pressure shift/Av 
isotope Buffer gas (cycles/sec mm Hg) (mm Hg) x 10° 

H! A -3.59 -2.53 

H! Ne +1.62 +1.14 

H! H,! -0.24 -0.17 

H! He +2.83 +2.00 

H’ A -0.63 -1.94 

H’ Ne +0.44 +1.34 

H3 A -1.98 -1.31 
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Table Il. The measured hyperfine splitting obtained by extrapolating to zero pressure. 








Number of bulbs Number of Hyperfine 

Hydrogen Buffer at different observations splitting 

isotope gas pressure per bulb in Mc/sec 
H! A 4 40 1420. 405712 
H! Ne ~ 20 754 
H! H, 4 40 721 
H! He 4 10 698 
Average 1420. 405726 
H? A 3 20 327. 384352 
H? Ne 4 15 347 
Average 327.384349 
H® A 3 40 1516. 701382 
H? Ne 3 40 410 
Average 1516. 701396 





the lines determined for the various buffer gases 
fail to meet in a point at zero pressure has been 
taken as a measure of the error in the hyperfine 
splittings of these two isotopes. The results for 
hydrogen and tritium are 


Av (H’) = 1420.405726 + 0.000030 Mc/sec, 
Av (HS) = 1516.'701396 + 0.000030 Mc/sec. 


The result for hydrogen agrees very well with 
the most recent atomic beam result of Kusch* 
and the paramagnetic resonance result of Wittke 
and Dicke® with no pressure shift correction. 
The value for deuterium is slightly outside the 
quoted error in the atomic beam measurement 
of Kusch.* The value for tritium disagrees with 


the atomic beam measurement of Prodell and 
Kusch.® 
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NUCLEAR RESONANCE ABSORPTION AND NUCLEAR ZEEMAN EFFECT IN Fe®”? 


G. DePasquali, H. Frauenfelder, S. Margulies, and R. N. Peacock 
University of Illinois, Urbana, Illinois 
(Received December 14, 1959) 


Mossbauer recently demonstrated that the 129- 
kev gamma ray in Ir’** underwent resonance ab- 
sorption when the osmium source and iridium 
absorber were both at low temperatures.’ The 
resulting absorption peak displayed the natural 
line shape.? Mossbauer explained this effect by 
adapting a theory of Lamb.* If the emitting atom 
is bound, at a low ambient temperature, ina 
solid of high Debye temperature, it cannot al- 


ways recoil freely. In a certain fraction of the 
transitions, the recoil momentum is given to the 
entire solid. The emitted photon then suffers no 
Doppler shift, but possesses the full transition 
energy and can be resonantly absorbed or scat- 
tered by another nucleus of the same type, also 
bound in a solid. Mossbauer’s experiments have 
been repeated and extended*»* and the theory has 
been refined.®>” 
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In the present Letter, we report on a resonance 
absorption experiment with Fe*’ in ferromagnetic 
iron. The pertinent nuclear data appear in 
Fig. 1(a).° The low recoil energy due to the 14.4- 
kev photon, 0.002 ev, and the high Debye tem- 
perature of iron, 335°K, permit experiments up 
to temperatures of about 200°C. The large cross 
section at resonance,® o,~1.5x10"** cm’, renders 
experiments feasible with natural iron even 
though the natural abundance, a, for Fe” is only 
2.17%. Furthermore, the Zeeman splitting of 
the nuclear levels due to the hyperfine coupling 
of the nucleus with the electrons is larger than 
the natural line width f (f =h/7=4.6x10~° ev) 
and can be observed directly. 

Our experimental arrangement is shown in 
Fig. 1(b). The source consists of carrier-free 
Co*” in less than 1 mg of iron. The absorber is 
an iron foil, sandwiched between two Lucite 
disks. The absorber disk can be rotated over a 
continuous range of constant speeds. The com- 
ponent of disk velocity in the direction of the 
gamma ray thus yields a continuous relative mo- 
tion between source and absorber and provides 
for the Doppler shift necessary to sweep over the 
absorption curve. A collimating slit is cut in 
such a way that the velocity component varies by 
about + 10% over the area of the slit. A Philips 
beryllium-window thin NalI(T1) scintillation coun- 
ter is used to measure the 14.4-kev gamma ray. 
The resonance absorption is expressed by a 
ratio 


€(v) = [N(~) - N(v) ]/N(~). 


Here, v denotes the average component of velo- 
city of the absorber in the direction of the in- 
coming gamma rays, N(v) is the corresponding 





FIG. 1. (a) Decay scheme of Fe®’. (b) Experimental 
arrangement. The velocity component of the rotating 
absorber disk in the direction of the gamma rays pro- 
vides a continuous and adjustable relative motion be- 
tween source and absorber. 
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counting rate, and N(«) is the counting rate at a 
velocity which is so large that no resonance ab- 
sorption occurs (v=1 cm/sec). 

The maximum resonance absorption ¢€(0) as a 
function of temperature of the Fe*®’ source was 
measured with a steel absorber of thickness 
16.5 mg/cm? at an angle 6=30°. The curve «(0) 
stays essentially constant from 25°C to 200°C, 
then begins to decrease slowly, and approaches ‘ 
zero asymptotically at about 550°C. 

Figure 2 displays «(0) as a function of absorber 
thickness. The absorption curve is analyzed by 
a method used by the Los Alamos group.*»*° If 
f is the fraction of the 14.4-kev line which is 
resonant, and f’ the probability of absorption i 
without recoil, then 


€(0) =f[1 -1,(x/2) exp(-x/2)]. 


Here, J,(z) =J,(iz) is the Bessel function of order 

0 and imaginary argument, x=/’nao,, and n is the 

number of iron atoms per cm’. Fitting the ex- 

perimental points to curves ¢«(0) vs x, calculated 

for several values of f, yields unique values for 

f and f’. From Fig. 3 it follows that f =0.20+0.04 

and f’=0.22 4 0.06. 
Figure 3 exhibits the resonance absorption 

curve e(v), not corrected for absorber thickness 

or finite velocity resolution. To discuss Fig. 3, 

we assume initially that the magnetic moment 

i, of the Fe*” ground state is zero and that the 

excited state in both source and absorber is split 

into four equally spaced magnetic sublevels. 

This splitting is due to the interaction of the mag- 

netic moment u* of the excited state with the 
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FIG. 2. Resonance absorption €(0) as a function of 
absorber thickness. The abscissa x =/f’ndo, is pro- 
portional to the absorber thickness. 
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FIG. 3. Resonance absorption curve e€(v) of the 
14.4-kev gamma transition in Fe*’, exhibiting Zee- 
man splitting. Source and absorber consist of natural 
iron and are at room temperature. 
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magnetic field H at the nucleus. The 14.4-kev 
lines emitted by the source then consist of four 
components, separated in energy by »*H/I*, 
where J* =3/2 is the spin of the excited state. At 
zero relative velocity between source and ab- 
sorber, these emission components match the 
corresponding four absorption components yield- 
ing the resonance absorption €(0). Moving the 
absorber with a velocity component v shifts the 
absorption components with respect to the emit- 
ted ones by an energy (v/c)E, where E =14.4 

kev. One can thus trace out the central line. 

Since both emission and absorption components 

are Breit-Wigner curves of width I, the central 
line should approximate a Breit-Wigner curve of 
width 2F=1.1x10-% ev. At velocities v=t+cyu*H/I*E, 
three of the four components should overlap again 
and produce, on either side of the center, a line 

of smaller intensity. Similarly, two more lines 

of even weaker intensity should appear on either 
side at twice and three times the distance of the 
first one. The ratio of intensities of the lines 
depends on the relative orientation of emitting 

and absorbing nuclei. For random orientation, 

one expects intensities 4:3:2:1. 

Figure 3 clearly shows the Zeeman splitting of 
the 14.4-kev transition. In order to interpret 
Fig. 3, some information about H is required. 
No direct measurement is available, but from 
the isotropic hyperfine splitting of iron group 
salts, a value H~=10° gauss is expected.’ 
This value will be used for the following discus - 


sion, although the actual field may be consider- 
ably larger.'* The assumption that the resolved 
lines are due to splitting of the excited state may 
be justified by considering the positions of the 
first lines that would result from a ground-state 
moment. With the experimental upper limit® 
Ug 0.05u,, and H = 10° gauss, these lines would 
appear at the positions of the arrows in Fig. 3. 
This ground-state splitting could not be resolved 
in our equipment, but it would broaden the ob- 
served lines if present. 

The measured value for the width of the central 
line, 2P,, =3.4x10°* ev, is considerably larger 
than the predicted one. This broadening may be 
due to a small magnetic moment of the ground 
state or to extranuclear interactions." 

From the separation of the observed lines the 
magnetic moment of the excited state may be 
estimated to be u*=0.50yu,,. Because of the un- 
certainty in H this value should be accepted with 
caution. A direct splitting of the lines in an ex- 
ternal magnetic field will yield a more reliable 
value of p*. 

During the present experiment, we enjoyed 
stimulating discussions with P. Axel, J. Hether- 
ington, J. D. Jackson, H. J. Lipkin, C. P. 
Slichter, and J. C. Wheatley. We would like to 
thank R. N. Lee for preparation of the absorber 
foils and J. B. Reynolds of Washington Univer- 
sity, St. Louis, for the production of the Co*”’. 
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POLARIZATION OF NUCLEAR RESONANCE RADIATION IN FERROMAGNETIC Fe*’t 


G. J. Perlow, S. S. Hanna, M. Hamermesh, C. Littlejohn, D. H. Vincent, 
R. S. Preston, and J. Heberle 
Argonne National Laboratory, Lemont, Illinois 
(Received December 28, 1959) 


The ferromagnetic property of iron has made 
possible a simple measurement of the polarization 
of resonance radiation from the level in Fe® at 
14 kev and a direct observation of a strong cor- 
relation in direction between the magnetization 
in a ferromagnetic domain and the internal mag- 
netic field at the nucleus. Recent observations 
have shown!‘ that the 14-kev radiation undergoes 
a strong absorption of the Mossbauer type’ at 
room temperature and that the absorption spec- 
trum displays? a hyperfine structure correspond- 
ing to the splitting of the 14-kev level by the in- 
ternal field. The present experiment consists in 
observing the resonant absorption with a ferro- 
magnetic source and absorber magnetized par- 
allel or perpendicular to each other. 

The apparatus is shown schematically in Fig. 1. 
On one side of the alnico magnet M is fastened 
the yoke structure designated by Y. The pole 
pieces P bridge the circular gap formed by the 
yoke, and the source S is mounted in the gap 
between the pole pieces. An identical structure 
for providing a magnetic field in the absorber A 
is mounted on the other side of the same magnet 
M. The design allows the pole pieces P and P’ 
(shown also in side view in Fig. 1) to be rotated 
relative to each other so that the relative orien- 














FIG. 1. 
resonance radiation in Fe*’, 
described in the text. 
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tations but not the magnitudes of the fields in the 
source and the absorber can be changed. In this 
arrangement the magnetic fields are perpendicular 
to the direction of propagation of the radiation 
which is detected with a thin (20-mil) Nal crystal. 
The source was in the form of an alloy of active 
Co” and Fe (99.91%) electroplated upon a strip 
of ordinary iron, 0.5 mil thick. This iron strip, 
attached to a piece of Lucite for rigidity, then 
bridged the gap between the pole pieces (see the 
side view in Fig. 1). The absorbing foil, a piece 
of rolled Fe*” (76%) metal, 0.1 mil thick, was 
mounted in similar fashion to the second set of 
pole pieces. The magnetic fields in the two gaps 
between the pole pieces were approximately 800 
oersteds, a value which should produce fairly 
good saturation in the foils. The photomultiplier 
of the detector was shielded magnetically and it 
was standard practice in a series of observations 
to rotate only the pole pieces P which are further 
from the detector. Under these conditions the 
influence of the changing field on the detector 
was observed to be negligible. 

The polarization of the resonance radiation was 
demonstrated by measuring the transmission 
through the absorbing foil as a function of the 
relative orientation of the fields in the source 
and the absorber. The result of one series of 
measurements is depicted in Fig. 2. It is seen 
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FIG. 2. Transmission curve of resonant radiation 
as a function of the angle between the magnetizations 
of the source and the absorber. The transmission 
obtained with a vibrating source and that obtained with 
zero field in source and absorber are also indicated. 
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that the transmission is strongly dependent on 
the angle between the magnetizations of source 
and absorber. Also indicated in Fig. 2 is the 
transmission obtained with a vibrating source 
and that observed when the fields in the source 
and the absorber were approximately zero. The 
latter measurement was obtained by replacing 
the magnet M with a nonferromagnetic substitute. 
The transmission with crossed fields is only a 
little less than that obtained with a vibrating 
source. On the other hand, with parallel fields 
the transmission is markedly lower than the “no 
field” value. The small deviation from symmetry 
about 0°(or 90°) that is apparent in Fig. 2 can 
possibly be traced to an instrumental asymmetry 
in the fields. As a check on the observations 
they were repeated with the fields in the source 
and the absorber reduced approximately to zero 
(see above). No significant angular effect was 
observed. 

The result in Fig. 2 is simply explained as 
follows. The alignment of the ferromagnetic 
domains in the source produces a corresponding 
alignment of the internal fields at the iron nuclei. 
Relative to this aligned field the radiation emitted 
in the transition from one magnetic sublevel to 
another is linearly polarized when viewed per- 
pendicular to the field. Hence the magnetization 
of the source serves as a polarizer of radiation. 
In like manner the magnetization of the absorber 
provides an analyzer of the polarization. 

This explanation has some interesting conse- 
quences. It provides a direct demonstration of 
the strong correlation in direction between the 
magnetization and the internal field at the nucle- 
us in a ferromagnet. It does not, however, give 


the sign of this correlation. In addition, it can 
provide information on the hyperfine structure 
of the levels in Fe®”. It has been suggested?» * 
that the splitting of the ground level, j=4, is 
relatively small (owing to the small magnetic 
moment of the state) and perhaps unresolved in 
the observed spectrum. A straightforward anal- 
ysis of the polarizations in the hyperfine spec- 
trum shows that the strength of the absorption 
in the case of a magnetized source and absorber 
will depend on the degree to which the ground- 
state structure is resolved, as well as on the 
character of the radiation (M1,E2). A quanti- 
tative interpretation of the data is complicated 
by the effect of absorber thickness and must 
await more detailed observations. 

Our thanks are due to D. R. Inglis andG. R. 
Ringo for interesting discussions; to M. R. 
Perlow for aid in the chemical preparations; to 
W. J. Ramler and A. A. Schulke for the cyclotron 
bombardments which produced the active Co*’; 
to F. J. Karasek for skillful preparation of the 
thin rolled iron foils; and to W. W. Managan for 
the preparation of the detector. 
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NEW “BREAK” IN THE O'%(y,)O** REACTION WITH A NEW SYSTEM OF MEASUREMENT 


D. Sadeh* 
Laboratoire de Synthése Atomique, Centre National de la Recherche Scientifique, Ivry, France 
(Received December 7, 1959) 


Using our 31-Mev B.B.C. betatron, we bom- 
barded cylinders of boric acid (H,B,O,) and 
measured the 2- minute g* activity of the O** with 
a new system. 

The cylinders (42 mm in diameter and 47 mm 
in length) were placed 25 cm from the betatron 
source of radiation and were irradiated for 5 
minutes; 20 seconds after the end of irradiation 
they were placed in a lead container between 2 





large Nal(T1) crystals (5 inches in diameter and 
4 inches in length). The two annihilation 7 rays 
were thus measured with a coincidence and 
single-channel analyzer system. In this way we 
used targets which were ~8 times thicker than 
those which are used when the BY is counted; the 
counting efficiency and solid angle are largely 
improved in comparison to previous systems 
where the 6* activity was measured by Geiger 
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reaction O'%(y, no". 


counters (~100 times more efficient than those 
of Penfold and Spicer,’ Bendel et al. ,? and Katz 
et al.*). The relative intensity was measured by 
using the same system to count the activity of 
cylinders of copper which were irradiated at the 
same time as the cylinders of boric acid. The 
activity of Cu® was counted for a 2-minute in- 
terval beginning 6 minutes after the end of the 
irradiation. As the activation curve of coprer 
is well known‘ in this region of energy, a correc- 
tion is made for each measurement. By meas- 
uring the intensity in this way and not with an 
electrical system which measures the current 
in an ionization chamber, we gained in two ways: 
(1) the error was 0.2% because of the large num- 
ber of counts; (2) fluctuations in the efficiency 
of the counters did not change the relation be- 
tween the counts from the oxygen and the counts 
from the copper cylinders. 

In this way we found all the “breaks” which 





were found by Katz°® and Penfold and Spicer, 
and another one placed between the threshold 
and the first “break” 110 kev above the threshold 
of the O'%(y,”)O** reaction (see Fig. 1). Cali- 
brating our betatron with the thresholds of the 
reactions F?*(y, m)F?® (10.40 Mev), Cu®™(y, 2)Cu® 
(10.78 Mev), O*%(y,)O**® (15.65 Mev),?*® and 
C(y,n)C* (18.72 Mev), this “break” is at 15.76 
+0.020 Mev. This is in agreement with the state 
that was observed in the N"°(p, ay)C™ reaction’ 
at 15.79+0.015 Mev. Thus, not only did we find 
a new excited state in the O'* nucleus but also 
this result indicates that calibrations® of beta- 
trons with the O'*(,,”)O** threshold were mis- 
taken by 110 kev (as was predicted recently by 
the authors of reference 6) by taking the first 
“break” as if it were the threshold. 

Our calibration was done by taking cylinders 
containing F’*® and Cu®, and measuring their 
activity in the way described above. This meth- 
od of calibration is so sensitive that even 5 kev 
above threshold we could see an increase of the 
yield. For the C(y,)C™ reaction we used a 
more powerful system which is described else- 
where.® Making the calibration in this way we 
cannot miss “breaks” which occur immediately 
above threshold, and thus we cannot make any 
mistake in determining the threshold. 
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PRODUCTION OF 2° AND A HYPERONS BY &* CAPTURE IN DEUTERIUM” 


Orin Dahl, Nahmin Horwitz, Donald Miller, and Joseph Murray 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received December 14, 1959) 


Recently, Ross reported a measurement of the 
branching ratio for the two reactions 


x +p—L°+n, (1) 
xD +p-A+n, (2) 


induced by the absorption of stopped 2 in liquid 
hydrogen.’ He found 


I(2°)/P(Z°+ A) = 0.33 + 0.05. (3) 


Day, Snow, and Sucher have tried to relate the 
analogous deuterium reactions 


x +d—D°+n+n, (4) 
Lt +d-A+n+n, (5) 


to those in hydrogen by means of an impulse- 
type calculation.” They pointed out that because 
of the small =~ - £° mass difference it is ex- 
pected that reaction (4) will be strongly sup- 
pressed. However, uncertainties in (a) the de- 
tails of the two-body = -p absorption process 
and (b) the strength of final-state interactions 
prevented a reliable estimate of the effect. 

An additional exposure of the 15-in. Alvarez 
deuterium bubble chamber to the separated K— 


beam of the Bevatron has resulted in an improved 
determination of the branching ratio for reactions 
(4) and (5).° A total of 145 charged A decays pro- 
duced by absorption of = ’s from the reaction 


K°+d-2 +nt+n (6) 


have been observed. The presence of the z* at 
the K” absorption vertex unambiguously indicates 
the production of a 2 , whether or not the range 
of the 2 was great enough to produce a visible 
track. 

In order to determine the A momenta, the co- 
ordinates of each event were punched onto IBM 
cards by means of the LRL precision measuring 
machine (Franckenstein). The events were least- 
squares fitted by suitable IBM programs which 
applied the energy-momentum constraints at the 
decay vertex. The accuracy of the analysis was 
checked by processing in a similar manner the 
monochromatic A’s from reaction (2). In Fig. 1 
we have plotted the idiogram calculated from the 
measured A momenta and their uncertainties. 

A recent remeasurement of the = -L° mass 
difference has been reported by Rosenfeld.* He 


found that M,.- -M,o= 4.45 0.37 Mev. Since the 
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FIG. 1. Comparison of observed spectrum with the prediction of the impulse model. 
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Q value for reaction (4) is then only 0.92+0.37 
Mev, the A’s are confined kinematically to the 


narrow momentum interval from 40 to 109 Mev/c. 


However, processes other than (4) may contri- 
bute to the observation of A’s in this momentum 
interval. We have considered the following 
sources of background: 

A. Spectroscopic analysis of the deuterium in 
the bubble chamber during the course of the ex- 
posure indicated that the molar contamination of 
hydrogen was ~2%. Ross found that A’s result- 
ing from decay of =°s produced in reaction (1) 
occurred in the momentum interval 0 to 135 
Mev/c. Therefore, a £° production in the resid- 
ual hydrogen in the present experiment would 
probably be indistinguishable from reaction (4). 
If we assume that = capture occurs in hydrogen 
in proportion to the amount present, one £° event 
and two A’s (of unique momentum 288 Mev/c) 
must be subtracted out. 

B. From Ross’s data it is also estimated that 
about five of the £ -d interactions occurred in 
flight. Because of the distribution in momentum 
of the = ’s from reaction (6), the in-flight events 
cannot be recognized. Since £° production pro- 
bably increases rapidly to the hydrogen value 
with increasing © momentum, we take account 
of the in-flight interactions in the calculation of 
the branching ratio by an additional subtraction 
of one £° and four A’s. 

C. A’s from reaction (5) may be produced in 
the expected momentum interval. Calculations 
were carried out in which the effects of possible 
final-state interactions were estimated. The re- 
sults suggest that the contamination of A’s from 
reaction (5) is negligible. 

We conclude that the branching ratio for reac- 
tions (4) and (5) when = hyperons are captured 
at rest in deuterium is 


PY, (2°)/T, (29+ A) = (5+ 3)/137 =0.0370.022. (7) 


This ratio is within the limits predicted by Day, 
Snow, and Sucher.? 

We estimate that the probability is less than 
10°* that seven or more events in the expected 
=° momentum interval could arise from the 
sources of contamination listed above. There- 
fore, the present result may be considered as 
providing an absolute lower limit to the = - £° 


mass difference: 


M,,- -M,o>3-53 Mev. 


It is of interest to determine the extent to which 


an impulse-type calculation is adequate to pre- 
dict the observed A spectrum. We performed the 
calculation under the following assumptions’: 

(a) =~ capture at rest proceeds from atomic S$ 
orbitals®; (b) absorption in the relative S wave of 
the = -p system dominates; (c) the effects of 
final-state interactions may be neglected. The 
Hulthén form was used for the deuteron wave 
function; the final A-2” system was described 
by a product of plane waves. Since the spin de- 
pendence of the two-body absorption operator is 
unknown, no attempt was made to include the ef- 
fects of the Pauli principle. The predicted spec- 
trum was found to be insensitive to the choice of 
relative 2-A parity. The result of the calcula- 
tion is shown in Fig. 1, with the known momen- 
tum resolution folded in. (In this comparison it 
is assumed that the effect of the four A’s from 
in-flight =~ captures is small.) The agreement 
is remarkably good. 

We conclude that both the branching ratio and 
the spectrum are qualitatively understandable 
in terms of a simple impulse model. There is 
no evidence (within experimental error) for dis- 
tortions of the spectrum due to the presence of 
strong final-state interactions. 

We wish to thank Professor Luis W. Alvarez 
for his advice and encouragement, and members 
of the bubble chamber and scanning staff for 
their assistance. 
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TESTS OF CHARGE INDEPENDENCE OF NUCLEON-NUCLEON INTERACTIONS* 


G. Breit, M. H. Hull, Jr., K. Lassila, and K. D. Pyatt, Ir 
‘ich Yale University, New Haven, Connecticut 
(Received December 14, 1959) 


the Evidence regarding near equality of specifically p-p and -p searches, respectively. 
nuclear p-p and p-n interactions arose from the The first test is similar to that quoted’? but 
consideration of ‘S states. Recent information’? differs from it in that same procedure rather 
of | regarding near equality of interactions is pro- than different procedures was used for n-p and 


vided for higher L in the comparison of values of p-p data and through inclusion of most available 
the Ps-ps coupling constant g,”? derived from the data. 


one-pion exchange potential (OPEP) adjusted to Some of the values of the unrenormalized 
p-p and n-p scattering data. Two tests concerned coupling constant are as in Table I. 
: with L >0 states have recently been performed in Error estimates without researching assume 
5s ' connection with an over-all phase shift analysis*® correctness of non-OPEP. The effect of re- 
> for scattering data with bombarding energies searching is seen in row 5 and is moderate. 
c- 9 Mev<E<345 Mev. It starts with phase shifts Rows 2 and 3 show stability to choice of OPEP 
of corresponding to approximate potentials‘ and re- parameters included in g,? search and increase 
E sults® of phenomenologic analysis at 310 Mev. It in normal accuracy resulting from extension to 
P improves over-all fits by gradient digital machine lower L andJ. Rows 6,7 and 9,10 show stability 
t : searches in phase-parameter space by factors of to independent variation in T=0 and 1 groups for 
~10 to 20 in D, the weighted sum of squares of n-p. Comparison of rows 6, 7, 8 with 9 and 10 
t deviations divided by the number of observations. and of 1 with 2 and 3 shows stability to refine - 
Roughly 400 and 300 measurements are used for ments in treatment of OPEP. Comparison of 


Table I. Values of the unrenormalized coupling constant Ly? from P-p and n-p data. 





Values used in search 





Run® Nucleons L and J By? Remarks? Row 
" YLAM p-p L>=5 1 13.5+0.9 a 1 
YLAM p-p J>3 1 14.0+0.2 b 2 
YLAM p-p L=5 1 13.3+0.7 b 3 
. YLA p-p L=5 1 12.5+0.8 b 4 
YLA p-p L>=5 1 12.0+1.1 b,R 5 
YLAN2M n-p L=5 0,1 14.5+0.7 a 6 
, YLAN2M n-p L>=5 0 14.541.1 a 7 
3 YLAN1 n-p L=5 0,1 14.1+0.8 a - 
YLAN1 n-p L>=5 0,1 13.2+0.8 b 9 
YLAN1 n-p L=5 0 13.0+1.0 b 10 





aKey to run designations: YLA as in reference 4; YLAM obtained from YLA by replacing phase parameters for 
*F;, °F,, coupling parameter to °H,, 'G, (the 3-4 group) by OPEP values, type 5, and releasing these for final 


to ; search; for YLAM(a) the 3-4 group was used with values from YLAM(b); YLAN1 obtained by gradient search start- 
ing with values of Gammel, Christian, and Thaler‘ with special treatment of K. 3 (first included in general search, 
. F then returned to the Gammel-Christian-Thaler value and other parameters researched, then K; released for final 


search); YLAN2M obtained by adjusting ¢ in gradient search with 6(YLAN2M) = 5(YLAN1) +(1-£)6(YLAN2) which 
gave §=1/4; YLAN2 uses in search start unpublished values privately communicated by John A. Gammel; in all 
n-P runs T=1 parameters are from YLAM. 

bkey to “ remarks’: @ ... relativistic one-pion treatment, accurate pion masses; 5 ... nonrelativistic one-pion 
treatment, pion masses approximated by m(n*); R ... non-OPEP researched in By" variation. 
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rows 6, 7 with 8 shows stability of g,* to varia- 
tion of search procedure for non-OPEP para- 
meters. Comparison of rows 1 to 5 with 6 to 10 
combined with that of 6 with 7 indicates charge 
independence of the OPEP. The mean of “a” 

type n-p values (rows 6, 8) is 14.3+0.5 and the 
mean of rows 1, 6, 8 is 14.0+0.5 while f? =0.08 
gives g,’ =14.4, 15.4, respectively, for m(m*) and 
m(n°). 

For the second test the m-p non-OPEP phase 
parameters corresponding to T =1 were varied 
from their p-p values which were initially used 
in gradient searches for T=0 parameters. In 
the first variation of the second test, D/D 
was calculated for T=1 parameters* 
b=Syp a+ E(OvrA - 5yp) with results as in 
Fig. 1. The difference in the best — is seen 
to be very small. The phase parameters 
varied were all those with L <5. In the second 
variation the 6 for T=1 were changed along 
the gradient of D according to 5=6°+ é’V,6, 
where v,5 is the gradient in phase-parameter 
subspace for T=1. Similar variations of T=0 
parameters were made employing direction 
V,0- Three calculations were carried out with 
search correction functions correlating changes 


min 








p-p case 








- 
2 





FIG. 1. Comparison of variations of D/D nin for 
n-p andp-p. This quantity is plotted against the para- 
meter ~ such that a change of — by 1 represents the 
difference between end points of searches employing 
different starting sets. 
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in the 6 for different E picked to have maxima 
at energies E = 50, 135, 300 Mev emphasizing 
different energy regions (“pk” for “peak”). 

For £ = 50 Mev the points for T=0 and T=1 

fall closely on the same graph, thus agreeing 

with charge independence, and D,,jn occurs 

for nearly the same phase parameters as for 

b-p. The probable error of the n-p determina- 

tion of phase parameters determined by the 

error matrix is high and the graph is not re- | 
produced. In Fig. 2 there are shown the T=1, 
n-p values of D plotted against &’ for E ox = 135 
and 300 Mev. Displacements of minimum from 
€’=0 are 0.0035, 0.0045, respectively, and 
represent the fractions 0.22, 0.11 of changes 

in &’ required for a change® of D by 1. Standard 
deviations S in phase parameters computed by 
the statistical matrix and errors C correspond- 
ing to AD =1 in the gradient cut give (log,,(S/C)) 
=7.1 and 3.6, respectively, at Ek =135 and 300 
Mev. In this somewhat arbitrary sense the dif- 
ference in p-p and n-p, T =1 phase parameters is : 
0.22/7.1 =0.11/3.6=0.031 of the standard de- 
viation of the determination of the parameters. 
The number of phase parameters was 11 and 
each contained 3 adjustable fitting parameters. 
The large error limits obtained by the statis - i 
tical matrix are caused partly by the large 
number (33) of fitting parameters. In fact the | 
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FIG. 2. Variations of D with distance along the 


gradient cut in T=1 subspace for Ex 135 and 300 
Mey in the ”-P case. 
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weighted sum of squares of uncertainties of 

the reproduction of measured quantities by the 
fit as calculated by the statistical matrix may 

be shown to be proportional to the number of 
fitting parameters. If the fitting were in terms 
of a theory containing only a few parameters, 

the fitting uncertainty would therefore be de- 
creased. For one adjustable theoretical para- 
meter producing variations along the search 
gradient used in Fig. 2, the standard deviations 
for E,,,=300 Mev computed from the variation 
of D with &’ would be 4,6=0.00033, 0.00054, 
0.0017, 0.0045, 0.0013, 0.0056 radian for 'S,, 
’P,, °P,, °*P2, ‘D2, °F, respectively, with sim- 
ilar small errors of remaining phase parameters 
for L<5. The largest uncertainty is for *F,. 
Similar small uncertainties apply to Eox= 135 
Mev. Differences between n-p and p-p phase 
parameters corresponding to the displacements 
At’ =0.005 and 0.0035 of D,,jy from their p-p 
positions are 1.8 and 4.1 times the uncertainties 
at E,),=300 Mev and 135 Mev, respectively. It 
it thus possible that charge independence is only 
approximate, with violations of the order of 2 
times the uncertainties in phase shifts listed. 
There is no assurance that the values of the phase 
parameters from gradient searches are as accu- 
rate as the A,6 suggest and therefore no definite 
indication of violation of charge independence is 
implied. The conjecture that the phase para- 
meters have been determined more accurately 
than the statistical matrix indicates, and that 
they may be of the order of the 4,5, appears 
more natural if it is remembered that gradient 
directions used for Fig. 2 were the last of a 
series of searches along some 12 previous gra- 
dient directions. They correspond, therefore, 

to directions of relatively gentle variation of D. 
It is improbable that a fundamental theory para- 
meter would give phase-parameter variations 
along the first few steep gradients but the chance 
of its giving variations similar to later gradients 
is greater. It is believed, therefore, that the 
absence of a large difference in n-p and p-p 








interactions suggested by Figs. 1,2 is not fortu- 
itous. 
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GREEN’S FUNCTIONS IN THE THEORY OF MANY-FERMION SYSTEMS* 


M. Bolsterli 
School of Physics, University of Minnesota, Minneapolis, Minnesota 
(Received December 7, 1959) 


The application of field-theoretic methods to 
the many-body problem’ has led to the use of 
Green’s functions in analyzing the lowest excita- 
tion energies of the many-Fermion system. In 
this note it will be pointed out that the usual 
single-particle Green’s function does not give 
direct information about the lowest excitation 
energies if there is a pairing energy of the type 
familiar in nuclear physics. An alternate Green’s 
function which will provide such information is 
suggested. 

The usual single-particle Green’s function is 


G, t; &’, t’) =i( THR, thy TR’, t”))), (1) 


where the expectation value is taken in the physi- 
cal vacuum, i.e., the lowest N-Fermion eigen- 
state of the total Hamiltonian. Here y and yt 

are the Heisenberg operators given by 


ue, =e ty Ge", (2) 


» &) =2),2,¢,(%), (3) 
where a, is the annihilation operator for a 
Fermion in the single-particle state e,(x), and 
the set e p(X) is assumed to be complete. Sub- 
stitution of (2) and (3) into (1) and use of the 
definition of the 7'-product gives 


G&, %’,7)= De Re, *(R)G(k,k’, 7), (4) 
k,k’ k k 


G(k, k’, T) 
aa i(E, -& 5)T ve t 
=i (Ola, |) la, 10), 7>0 
-i(E, -Ej)t 


=D ¢ (Ola, ,"12)<tla, 10), 7<0, 


(5) 


where 7 =/-?’, the sum j runs over all states of 

the (V+1)-Fermion system, and the sum 7 runs 

over all states of the (N - 1) - Fermion system. 
In general, k is a conserved quantity (k is 


momentum or angular momentum), so that 
a tT \ 12 
| = | js 

(0 a, |j)<jla,, 10) 6, pr! COle, >| (6) 


Then 


G(k, k’, 7) =6 


k, poles T), (7) 


and G(k, 7) can be written in the form 





G(k, r) = Limca) -f 7 G(k, a, ne "an, (8) 
with? 
(Ola, 17) I? Ola. Tz) 1? 
G(k, d, ae? . (9) 
A+E, ? grin SN - Eg +L) -i0" 


Now suppose that there is a pairing energy 


term in the ground-state energy of the N-Fermion ° 


system, i.e., that 

E,(N)=LN, N even 

E,(N) = nh +4,, N odd (10) 
where €p is independent of N for large N. It is 
clear that Ep is negligible insofar as it affects 
E,/N, the energy per particle, for large N, but 
we will see that it is not negligible insofar as it 
affects the Green’s function G(k, A, n). From (10) 
it follows that 


E,(N) ~B AN +1) =-p F ey ~€ (N+ 1), 
E,(N)-E (N -1) =u Fe, -€(N-1), (11) 


where the upper (lower) sign is for even (odd) 
N and €;, €; are the excitation energies of the 
states j,l, 


€(N+1)=5 (N+1)-E (N +1). (12) 


Substitution of (11) into (9) gives 
1¢0| j) |? 
(0 a, \j) | 





G(k, A, n=-D 
j A-H sie” 
(Ola, Fiz) \2 


Shs stdin 





- a (13) 
Fj A= PSE, +6 MN -2)-89 


where, again, the upper (lower) sign is for even 
(odd) N. If €p is positive, corresponding to the 
pairing energy observed in nuclei, then the 
zeros in the A plane of the denominators of 
G(k, 4, 1) for all k are shown in Fig. 1, where 
€j,€) are assumed to be continuous with €; min 
=¢€ min = 0. 

Since G(k, A, 7) is usually calculated for N even, 
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{ Fig. 1(a) shows that the existence of a gap be- pairing energy is less than half this gap. For 
- tween the particle and hole parts of the spectral noninteracting Fermions the pairing energy is 
representation of the single-particle Green’s easily seen to be zero. 
function does not necessarily mean that there It is also clear, however, that if one can go 
is a gap in the excitation spectrum of the sys- beyond the lowest excitation with wave number 
tem. The gap in the spectral representation may k and find the locations of the two or three low- 
be due to a pairing energy €», while €; may be est poles of G(k, A, n) for a fixed k, then a gap 
(8) continuous starting at «;=0. In order to prove in the excitation spectrum can be inferred from 
: the existence of a gap in the excitation spectrum, a gap between the first and second poles of 
it is necessary to show that the spectral repre- G(k, A, n) for k near the wave number for which 
sentation of the Green’s function has a gap be- the excitation energy is minimal. 
tween its particle and hole parts and that the In order to avoid the difficulty with the pairing 
“» (9) energy, consider the Green’s function® 
= , Keke » (ode 2. Hoe (+) 
| G% t; x’, t’) =i(T (a Me thy (x, t)y Mgr, t’)y gr, t’))). (14) 
mion | 
Here the superscripts + and - in parentheses indicate, respectively, the creation and annihilation parts 
of an operator with respect to the Fermi sea: 
10) 
' oo iHt -iHt - 
. 4 iG, =e > ae (Xe ai YG, =e! > ae et ; (15) 
| k<k 
3 F mk, 
ut 
it T > = - -j),> + > 
10 Oe, 9 <4 NG, 0, OG, =O, 0. (16) 
Now, application of (14) gives 
> — * > > * > - , 
11) | G%, t,t”) = 22h Be (Xe, (Xe, (x')e, (%')G (Rey, Ras gy Rast -t ), (17) 
ky, ky <R ps ka, kg? RE : . : 
a G (Rikakgka; T) 
€ 
12) — = jetHoT (a A o tht, . 908 
in k, “ky hy “hy 
“iE , 
(a) zie "(a Pe edd ', dig REM (18) 


ky hy k, hk, 


We now specialize the case where the ep(x) are 
13) _ plane waves V~“*exp(ik-x). Then the conserva- 
ie tion of momentum requires that 











en mn 
; k, -k, =k, -k, =q. (19) 
(b) It then follows that 


FIG. 1. The zeros in the complex A-plane near A =p siete ht 1 (x - xX’) 
of the denominators of G(k,A,n) for (a) N even and Go (x, t; x’, t’) -15 2/4 Gp (G,¢-t’), (20) 
(b) N odd. A=p is not a zero. q 












en, 
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where 
G4, 7) =D es sae aad (Olas lj) 1?, 7>0 


=D eo FP \(o1ag"9) 1, <0 
(21) 
A.= ) a> a> - 
q lkI< kp k k+q, (22) 
lk+ql>kp 
T T 
An = } a> a> - 
q lk I> kp k k-q. (23) 
lk -qi<kp 


Note that the states j in (21) are states of the 
N-Fermion system. Then, as before, 


G7) = lim(2n)"'f _G, Ga, ne" aa, (24) 





7-0 
*\ 2 | Tt . 2 
po fiOidg tdi? OiAg ip 
G (4, A; =F} h-€.+i — ioe. <i 
J j j 
(25) 


Thus except for states which are not created by 


—— 


Az or Ax' the spectral representation of G;, 
gives just the excitation energies, «€ j of the 
system, and it might, therefore, be called the 
excitation Green’s function. Although the numer. 
ators of the two terms in (25) look the same, 
they are not, since the restrictions on the k’s 
in (22) and (23) are different. In fact, it is easy 
to see that for noninteracting Fermions only the 
first term contributes. It is interesting to note 
also that €, and » can be obtained by comparing 
just the lowest poles in the spectral representa- 
tions of the two Green’s functions, G and G,. 





* 
Supported in part by the U. S. Atomic Energy Com- 
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34, 417, 433 (1958) [translation: Soviet Phys. JETP 
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5It is also possible to use the same function without 
the - and - superscripts; this is the procedure used 
by DuBois (reference 1) in analyzing plasmons. 





UNITARITY CONDITION BELOW PHYSICAL THRESHOLDS 
IN THE NORMAL AND ANOMALOUS CASES* 


S. Mandelstam 
Department of Physics, University of California, Berkeley, California 
(Received December 7, 1959) 


It is often the case with dispersion relations 
that the absorptive part for positive energies 
extends below the physical threshold. Thus, for 
example, in the process 21 ~N+N or »-N+N, 
the physical threshold is at t =4M? (¢ being the 
square of the energy), whereas the absorptive 
part extends to ¢=4,,", the square of the energy 
of the lowest intermediate state. When calcula- 
ting the amplitude for these processes,’ the uni- 
tarity condition has to be used in the region 
4u?<t<4M”’, and the question at once arises 
whether it is really valid there. Indeed, for 
those processes where there are anomalous 
thresholds, the imaginary part of the transition 
amplitude is known to be nonzero in a region 
where it should be zero according to the unitarity 
condition. It has therefore not been possible up 
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until now to apply dispersion relations to the 
calculation of such processes. The object of 
this Letter is to attempt to justify the use of the 
unitarity condition below the physical threshold, 
and to investigate how to modify it in the anom- 
alous case. We shall find that the same treat- 
ment is applicable to the normal and anomalous 
cases —it will confirm that the unitarity condition 
as used in reference 1 is correct, and will also 
show how to calculate the imaginary part of the 
transition amplitude in the anomalous case for 
the entire region in which it fails to vanish. We 
shall assume the validity of dispersion relations 
for partial waves throughout, and shall also 
approximate the unitarity condition in the physi- 
cal region by taking only two-particle inter- 
mediate states. 
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First let us quote the results which are ob- 
tained in the normal case if the unitarity condi- 
tion is used below the physical threshold. We 
consider the process 21 ~N +N, and neglect the 
spin of the nucleon, which merely complicates 
the algebra. The transition amplitude for the 
§ wave (any other angular-momentum wave could 
be similarly treated) may then be written as 


A(t) =A™(t) +A(2), (1) 
where - 
m ,ImAM(t’) 
amy =2f at par aa (2a) 
ImA(t’) (2b) 


1 oe 
A(t) -1{ tp 
T 4p? t -t 


The limit a is less than 4”. The absorptive 
part ImA” comes from pion-nucleon scattering 
and will be assumed known. It is then required 
to calculate A”, whose absorptive part comes 
from the reaction in question, by unitarity. 

In most of the calculations performed up till 
now, the perturbation result for A(t) is taken. 
The contribution of the nucleon poles to the 
scattering amplitude (not just the S waves) is 


s-M? u-M?’ 





(3) 


where s and uw are the squares of the energies of 


_ the pion-nucleon and the crossed pion-nucleon 


systems. Expressed in terms of the energy and 
angle for the reaction 21 —~N +N, 


s =2pqcosé -p” -¢, (4a) 
u = -2pqcosé - p? - q’, (4b) 
pand g being the momenta of the nucleon and 
meson, given by 
p® =3t-M’, (5a) 
q =3t-p’. (5b) 


On inserting these expressions into (3) and iso- 
lating the S wave, we find that 


wy). £0? -¢ -M? + 2pq 

OO 9g pF MF -2pq" ” 
As long as the signs of the square roots are 
defined to be the same inside and outside the 
logarithm, it is immaterial which sign is chosen. 
We take that branch of the logarithm which is 
real for ¢>4u?. As q is imaginary and p real 
for 4u?<¢<4M?, (6) is often written as an in- 
verse tangent. 





The unitarity condition for ImA® in the two- 
pion approximation is 


ImA®(M’, t) 
=(q/8rw){A™(M’, t)+A(M’, t)}B(t), (7). 


where B(t) is the S-wave amplitude for pion- 
pion scattering. The arguments M’ are inserted 
for future reference and should be ignored now. 
w is the center-of-mass energy vt. Equations (7) 
and (2b) can now be solved for A®. This has 
been done by Frazer and Fulco’ and by Omnés.? 
Either form is equally suitable for us at the 
moment, but the Omnés solution will be slightly 
more convenient in the subsequent discussion 
and we shall adopt it here. The result is then 


in ap LNNAM(M', t’) 
4p? w'(t’ -t) ' 








A®(M’, t) = (8) 


1 
87°D(t) 


D and N are defined by 


D(t) = exp} tf au 2 . 
4u? 


(9a) 


N(t) =D(t)B(t) 


98% sind exp pif at Ue) , (9b) 
q 7 4p! t'-t 


6 being the pion-pion phase shift. For ¢ real, 
D(t) is defined in the usual way by letting ¢ ap- 
proach the real axis from above. The phase 
factor of Dit) when ¢ is real and greater than 
4u? is e~*©, so that N(t) is real and has no cut 
over this range of values of ¢. 

It is now necessary to justify Eq. (8) without 
using the unitarity condition below the physical 
threshold. To do this we consider the Green’s 
function (in momentum space) rather than the 
scattering amplitude for our process. The 
Green’s function is defined for arbitrary values 
of the squares of the momenta of the particles; 
if we multiply it by i and by a factor (p? - m’) for 
each particle, where m is the relevant mass, 
we obtain a quantity which is equal to the scatter - 
ing amplitude when the squares of the momenta 
are equal to the squares of the corresponding 
masses. We shall, however, allow the squares 
of the nucleon momenta to have an arbitrary 
value M” (the same for both nucleons); we can 
then isolate the S wave and obtain a quantity 
A(M’,t), which is equal to A(t) when M’=M. 

By the usual technique of evaluating the ima- 
ginary part of the Green’s function, it is then 
easy to see that A(M’,t) satisfies the unitarity 
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equation (7) if we are in the physical region, i.e., 
if t>4M”. Instead of (6), A(M’, t) will be given 
by 


-p’? -  - M*+2p'g 


g*(M’) 
; inp” -@ -M -2p"q’ 


b’q 





AM’, t) = (10) 


where 


p” =3t -M’. (11) 


Note that M? in (6) is not replaced by M”, as it 
refers to the mass of a real intermediate state. 
g’(M’) is the square of the vertex function when 
the square of the momentum of one of the nu- 
cleons is M”. 

If now M’<yu, the whole range of values ¢> 4? 
where ImA®(M’, ¢) is nonzero is in the physical 
region, so that Eq. (7) is always applicable. We 
can therefore solve it and obtain Eq. (8) for 
A®(M’',t). This solution is thus justified if 
M’<u. 

To obtain the physically interesting case where 
M’=M, we now make an analytic continuation of 
this solution as a function of M’. This is per- 
missible, as the Green’s function in the physical 
region is known to be the boundary value of an 
analytic function.* We therefore have to examine 
the analytic properties of A™(M’,t), given by 
(10). The quantity 2*(M’) is analytic if M’<(M+ uy)’, 
since it satisfies a dispersion relation in M”, 
which can in fact be proved.* The remainder of 
the expression will have a branch point where the 
numerator in the logarithm vanishes, i.e., at 


t =t, = -M? +2(M” + yp?) -(M” - p?)?/M?, (12) 


This value of ¢ is always less than 4” unless 


M” =M?* + p?, (13) 


when ¢, =4u”. If M”<M?+ ", the right-hand 
side of Eq. (10) has no branch point at p’=0 or 
q=0. 

Equation (8) can now be analytically continued 
from values of M’ less than p to M’=M. The 
equality (13) is then never satisfied, so that the 
singularity at the point given by (12) is always 
at a value of ¢ less than 4p”, i.e., below the 
range of integration in (8). There is also never 
a singularity at g=0. The right-hand side of 
(8) is accordingly an analytic function of M” in 
the relevant range, so that it can be continued 
onto the mass shell. Equation (8) is thus cor- 
rect if M’=M, and we get the same result as if 
we had applied unitarity directly, with M’=M, 
below the physical threshold. 

We next examine the calculation when there is 
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an anomalous threshold. Suppose the theory 
contained, in addition to the meson and nucleon, © 
a baryon of mass B, where 


M? > B? +. (14) 


Equation (3) would then contain terms with M 
replaced by B, corresponding to an intermediate 
state of this baryon in pion-nucleon scattering, 
Let us now calculate the amplitude for the pro- 
cess 21 ~N +N, and use these new terms in the 
calculation of AY. Equation (10) is thus replaced | 
by 
ayy 7) OM) 2b? -G - B* + 2p'q 

AM D2" ig 8 pha - BP - 2p" 
The solution (8) can be derived as before if | 
M’<u and, for these values of M’, there is no 
anomalous threshold. Equation (12) for the 
singularity becomes 





(15) 


t=t, = -B?+2(M" +2) -(M? -p2)?/B%, (16) 
t, is less than 4u? except when 
M? =B? + y?. 


(17) 
If M> B* +, (15) will behave like -2ni/p’g at | 
q=0. 

The analytic continuation in M’ is now not 
straightforward, as the value (17) is reached 
before M’=M and, when M” has this value, 
t,=4u", so that the singularity in A"(M’, t) is 
at the end of the range of integration in (8). In 
order to avoid this, let us give M” a small posi- 
tive imaginary part. As M’ increases, the 
singularity at ¢=¢, moves as indicated by the 
dotted line in Fig. 1, ultimately reaching the 
point ¢,, given by 

t, = -B? + 2(M? + py?) - (M? - ?)?/B?. (18) 
When M” was less than B?+?, i.e., before the 


singularity went round the point ¢=47, the 
integral in (8) was along the real axis from 4,” 





_——— ee ~~ 


t=t, t=4u2, 
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FIG. 1. The singularities in the integrand, and the 
path of integration, for the anomalous case. 
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to «. When M”>B*+,?, therefore, we have to 
integrate along the solid curve shown. Note that 
the integrand in (8) is an analytic function of ¢ 

in the relevant region, so that we are permitted 
to deform the contour —we have already pointed 
out that N(t) has no cut along the positive real 
axis. 

We may now let the curves in Fig. 1 approach 
the real axis. The difference between the value 
of (15) on the two sides of the cut for t,<t<4y? 
is -2ni/p’q, so that (8) is replaced by 


1 "ia Nit’) 

(2) , So oe hh a 
AM" =a Dw), Tp wit 
0 
q'N(t)A(mM’, 


t’) 
weer-t) °°) 





1 J 
a 
817D(t) 4p? 


AM’, t), given by Eq. (8) if M’?<B* +? and by 
Eq. (19) if M’ > B?+, is thus an analytic func- 
tion of M” if t¢>4u*. It might be thought that 
there is a branch point when M” = B+”, but 
this is not the case, as we would have obtained 
exactly the same result had we given M” a nega- 
tive instead of a positive imaginary part. The 
curves in Fig. 1 would then be reflected across 
the real axis, but the sign change in the differ- 


ence of the logarithm on the two sides of the cut 
would be cancelled by the change of sign of q’ 
in the integrand of Eq. (8). 

When M” > B? +7, we have seen that A"(m’, f) 
behaves like -27i/p’q at q=0, so that ImA(m’, t’) 
is equal to -g*(M’)B(4u?)/167u(M” - p?)”” instead 
of zero at t=4y?. The cut in A® will extend to 
the anomalous threshold at ¢=/,. By putting 
M’=M in (19) we obtain the result of physical 
interest, with an anomalous threshold at ¢ =¢,. 
The method of applying the unitarity condition 
only above the physical threshold, but for vary- 
ing masses, and then continuing analytically in 
the masses, can thus handle both the normal 
and anomalous cases. The extra cut in the ano- 
malous case is a mathematical consequence of 
the analytic continuation, and does not appear to 
have any precise physical significance. 





“Supported by the U. S. Air Force Office of Scientific 
Research. 

‘Ww. R. Frazer and J. R. Fulco, Phys. Rev. Let- 
ters 2, 365 (1959), and to be published. 

*R. Omnes, Nuovo cimento 8, 316 (1958). 

3G. Kallen and A. Wightman, Kgl. Danske Videnskab. 
Selskab, Mat.-fys. Skrifter 1, No. 6 (1958). 

4A. M. Bincer (to be published). 





NEW TEST FOR A/=1/2 IN K* DECAY 


Steven Weinberg 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received December 28, 1959) 


We wish to suggest a practicable test of the 
4J=1/2 rule, based on a comparison of the pion- 
energy distributions in 7* and t*’ decay. At 
present, the only check of AJ=1 /2 in these proc- 
esses is the successful prediction of the 7’/7 
branching ratio.’ However, it is well known 
that the branching ratio tells us that AJ=5/2 and 
4J=7/2 are absent but tells us almost nothing 
about the possible presence of AJ =3/2 terms. 

The only symmetric three-pion states have /=1 
or /=3, and the other, nonsymmetric and hence 
inhibited states with J=1 or J=2 (which could be 
produced by a AJ =3/2 term) cannot interfere with 
the symmetric states in a measurement of decay 
rates. It is of course very important to learn 
whether the nonleptonic weak interactions involve 
a mixture of AJ=1/2 and AJ=3/2. In particular, 
ithas been noted that such a mixture would re- 


sult if these interactions arose from a folding of 
a 4IJ=1/2 strangeness-nonconserving current 
with the usual AJ=18-decay current.? From ex- 
perience with the 7’/7 ratio we see that a test 
for 4]=3/2 terms must depend on measurements 
of pion asymmetries of some sort. 

Suppose we let A,(7,7,7,) and A,/(7,T,T;) be 
the Lorentz-invariant amplitudes for K* decay 
into 1*1*n~ or 1°n°n* with kinetic energies T,, 
T,, and T,, respectively. The Bose statistics of 
pions implies that 

ATT oT) =A i 77. 


27173 (1) 


for j=7 or tT’. We shall break up A j into sym- 
metric and nonsymmetric parts: 


Ss N 
A (TT 973) =A, (7,T,7,)+4, (7,7,7,), (2) 


87 






VOLUME 4, NUMBER 2 


PHYSICAL REVIEW 





LETTERS JANUARY 15, 196 





S = 

A, (7,7,7,)= $A (7, 7,7.) +A (737 ,T,) 
N sis - 

A, (7,7,T,)= [2A (7,775) A(T 37,7) 


~A(T)T 37 )) (4) 


Now, the six amplitudes corresponding to 37 
states with definite total ] can be expressed as 
linear combinations of the six amplitudes ob- 
tained by permuting the arguments of A; and Az’. 
When we set the /=2 amplitudes equal to zero, 
we obtain 


P ae ee (5) 


Also, by setting the /=3 amplitude equal to zero, 
we have 


A_ =-2A_. (6) 


Hence we see from Eqs. (2), (5), and (6) that, 
under the AJ=1/2 rule, the energy dependence 
of the 7’ matrix element is entirely determined 


by the behavior of the matrix element for 7 decay. 


The presence of AJ =3/2 terms would invalidate 
Eq. (5), though not Eq. (6). 

To make this result more concrete, let us use 
the coordinates of the Dalitz-Fabri plot,* 


x =psing =V3(T7, -7,)/®, 
y =pcos¢ =(37, - Q)/&, (7) 


where we have Q=7,+T7,+T7,. We shall expand 
A in the series 


Q.\n 
=D a(n, a(t) p cosig, _— (8) 


n,l ins 


A(T.T.T 
J K 


1 2 3 


where a(n, /)=0 unless n -/ is even and non-nega- 
tive. [Equation (8) may be justified by performing 
a multipole expansion of A, which will yield a 
power series in |P,|?, |P,|?, |P,|. Since we 
have |P|?=2m7T+T?, this can be rewritten as a 
power series in T,, T,, and T,, and hence in 

Qx and Q@y. Since, according to Eq. (1), we can- 
not have terms odd in x, this power series can 
be written as the cosine series of Eq. (8). ] 

It may easily be seen that the terms in Eq. (8) 
with 7=0,3,6,9,... are completely symmetric 
and belong to AS, while the others belong to AN, 
Thus, according to Eqs. (5) and (6), we obtain 


a_(n,1)=a_{n, 1) (9) 


: 


for /=1,2,4,5,7,8,---, and 
a(n, l)=-2a_,(n, 1) (10) 


for 1=0,3,6,9,---. These equations are ex- 
pected to be correct up to Coulomb and mass- 
shift corrections of a few percent.® Since the 
a;(n, 1) are only slowly varying functions of Qj, 
the 12% difference between Q; and Q;, will also 
not alter Eqs. (9) and (10) by more than a few 
percent at most, though it will have a very large 
effect on the total available phase space. 

In principle, for each (n,/), Eq. (9) or (10) can 
serve as a test of the presence of AJ=3/2, 5/2 
terms or of AJ=5/2, 7/2 terms. There is one 
case, however, for which Eqs. (9) and (10) can : 
be tested by merely observing 7, in a reasonable : 
number of 7 and 7’ events, without having to 
measure any neutral pion energies. We expect 
theoretically that all of the a;(n, 1) should be 
roghly comparable in magnitude. [Very strong 
final-state interactions might give rise to some 
factors of (m,/m,)”", but this seems unlikely. | 
Since Qj/mK is small, we can try neglecting all 
terms in |A.|? except for |a;(0, 0)|? and the in- 
terference between a,(0, 0) and a,(1, 1). We then 
have 


Q? 
1A = + ss S \ 
A (7,775)! = af B{@./m dy + o(s)): (11) 





where 


@, = 1a (0, 0)1?, (12a) 
, 8 


and 


B ;=2Rela (1, 1)/a (0,0)}. (12b) 


If we apply Eqs. (9) and (10) to Eqs. (12) and (13), 
we obtain the familiar result @,,=@,/4, and the 
new prediction that 

B_,=-28_- (13) 
(It is of course helpful experimentally that |8,/! 
is twice |8,1|.) To calculate rates from Eq. (11), 
we must multiply by the relativistic phase-space 
density, which is conveniently a constant for 
both 7 and 7’ decays, and integrate over unob- 
served energies. We then obtain the differential 
decay probability 


= 14) 
w (y)dy [1+ BAQ,/m KY kk io )dy, ( 











where xj;(y), the maximum value of x for a given 
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y, is given exactly by 





(1+y)(C.+ny)(A,-By)7” 
x0)-| j_j__J | (15) 


1+37(A.-By) 
% J J 


Here we have Aj=1- 5j+Nj, Bj=1 + Oj +20; , and 
Cj= =1+35;+n;, where 6, =0, re =@,/6my, 
6, = (mq - mo)/3M a0; Nz =Q, 1/6m yo: 

It follows from Eq. (14) that the distribution 
in unlike and like pion energies will be 1+a;T, /mx 
and 1+y, .T,/mx, respectively [except for a phase- 
space gactor similar to Eq. (15)], and that we 
have 


3B. 
"TB /m,)° - 
eh 
and 
=-1 17 
¥;=~36,- (17) 


Ina recent study of 959 7 events,° it was ob- 
served that the deviations from “phase space” 
of the distributions in T, and in T, were roughly 
linear, with a,= 6.82 1.2, y; =-2.240.3, in fair 
agreement with Eq. (17) and hence with our neg- 
lect of quadratic terms in Eq. (11). We then have 
8;~1.3, in agreement with our expectation that 
the a,(n,1) should be of comparable magnitude. 
We predict 8,,~-2.6. The intrinsic uncertainty 
in this prediction, due to neglect of quadratic 
terms in Eq. (11), is probably about 20%. By 
careful study of existing 7 data, it should be 
possible to estimate a,(2,0) and a,(2, 2) [for 
example by checking (17) more accurately] and 
hence refine the accuracy of our prediction of 
the 7* energy distribution in 7’ decay.” 


It is a pleasure to thank Professors R. H. 
Dalitz and M. A. Ruderman for their comments, 
and Professor D. Glaser, Professor G. Gold- 
haber, and Professor S. Goldhaber for their 
suggestions on experimental possibilities. 





* 

This work was done under the auspices of the U. S. 
Atomic Energy Commission. 

'R. H. Dalitz, Proc. Phys. Soc. (London) A69, 527 
(1956). 

2Okubo, Marshak, Sudarshan, Teutsch, and Wein- 
berg, Phys. Rev. 112, 665 (1958). 

34 similar situation for K° decay has been discussed 
by S. B. Treiman and S. Weinberg, Phys. Rev. 116, 
239 (1959). 

‘R. H. Dalitz, Phil. Mag. 44, 1068 (1953); E. Fabri, 
Nuovo cimento 11, 479 (1954). 

‘This remark applies only to effects arising within 
the range of the strong interactions. In order to take 
account of the long-range Coulomb final-state inter- 
action in tT decay it is necessary to multiply the right- 
hand side of (8) by a Coulomb correction factor 
1+ (am/2)[(1/243) + (1/293) -(1/242)], a8 shown by Dalitz, 
reference 1. Empirical distributions should be divided 
by the square of this factor before analyzing to find the 
a(n, 1) coefficients. This was not done here, and the 
values quoted below for 8 are therefore too small. We 
wish to thank R. H. Dalitz, H. P. Noyes, and M. A. 
Ruderman for helpful discussions on this point. 

®McKenna, Natali, O’Connell, Tietge, and Varshneya, 
Nuovo cimento 10, 763 (1958). Of the events studied, 
419 were from Baldo-Ceolin, Bonetti, Greening, 
Limentani, Merlin, and Vanderhaega, Nuovo cimento 
6, 84 (1957). 

"R. H. Dalitz has performed an analysis of 900 
earlier t events [Reports on Progress in Physics (The 
Physical Society, London, 1957), Vol. 20; and private 
communication] and obtains 8, =1.6 + 0.5 and 
Re{[a7(2, 0) -@7(2, 2)]/a7(0, 0)}=-8.4 +6.5. 








HIGH-ENERGY LIMIT OF SCATTERING CROSS SECTIONS 


D. Amati, M. Fierz, and V. Glaser 
CERN, Geneva, Switzerland 
(Received December 14, 1959) 


Pomeranchuk’ has shown that under rather 
plausible assumptions concerning the very high 
energy dependence of total cross section—the 
main one being that they behave almost as con- 
stants at infinity—the difference of the cross 
sections for a particle and its charge conjugate 
on the same target tends to vanish at infinity. 
More explicitly: if o+(E) and o~(E) are the total 
cross sections for a particle and its charge con- 
jugate (r+ -7~, proton-antiproton, K* -K~, K°-K°, 


etc.) at total energy E of the incoming particle 
in the laboratory system on a specific target, 
and if 


lim o*(E)=0'(«), 


oO 


lim o (E)=0 (), 
E~o 


(1) 


o*(~) and o~(«) being finite constants (or zero), 
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then Pomeranchuk shows that 
0 («) =0 (~), (2) 


Pomeranchuk tries to justity (1) by pointing out 
that due to the finite range of interaction between 
particles, the quantity 


A* (E)/E, 


where A*(E) is the forward scattering amplitude 
for the corresponding process averaged over 
possible spins, should go to a limit with increas- 
ing E. Though it seems evident that the said 
quantity is bounded, it is not clear that if it does 
not go to zero it tends to a limit. But for Pom- 
eranchuk’s argument (1) it is essential that at 
least the imaginary part has a limit. We shall 
explore some consequences that can be drawn if 
the condition 


lim A’ (E) /E =constant (3) 
E — oO 


really is fulfilled. We note that by virtue of the 
optical theorem the condition (3) on the imaginary 
part of A+(£) is equivalent to (1); for the real 
part, D*(£), (3) implies 


: 
lim a =constant. (4) 
E-« 


The property (2) is of importance for the appli- 
cation of dispersion relations: in fact, in many 
cases where the unsubtracted dispersion rela- 





+ P . + . - 
D )=¥ (1455)? (M) +5 1-7 D ‘unsDz 


tions for the difference of particle and antipartic|: | 
amplitudes can converge because of the vanish. ~ 
ing of 0+ - 0” at very high energies. Such dis- 
persion relations, subject to that hope, were 
used and analyzed for several processes and 
proved to give meaningful results. This is the 
case for 7 -N S-wave scattering” (for the com- 
bination a, - a@,) and for K -N scattering.* In 
order that those applications make any sense, : 
it is necessary and sufficient, however, that 
integrals of the type 

to" (E) - o*(E) 

f ee (5) 
converge for t~«, Even if (2) makes more plau- 
sible the convergence of (5), it is by no means 
sufficient: to ensure it, it would be necessary 
to know how o~ - o* goes to zero at very high 
energies. 

We note that if 


lim [o (E) -o (E)]~(1/inz) (6) 
E-« 


(behavior which in some cases has been interpreted 

to be Pomeranchuk’s prediction), then (5) would 

diverge as InIn¢. In reference 1, however, no 

discussion is made on how the limit (2) is reached, 
We want to show that using the same starting 

point as Pomeranchuk, something can be said 

on how o~ - o* goes to 0; in fact we show in this 





note that the integral (5) is indeed convergent. 
Let us start from the subtracted dispersion 
relation whose convergence is guaranteed by (1): 


ot(E’) o7(E")\ dE’ 
af ‘ — Sra p’’ (7) 








where p =(E”+M?)”?, M being the mass of the particle in question and 2, ,App?/(E +E) the contribution 


of possible poles (bound states at energies Ey), Ap, being constants. The lowest limit of integration £, | 


' 


depends on the process in question (possible existence of unphysical regions). We could start from the 
dispersion relation for D~(E) without changing the conclusions we shall reach. 


Let us define 


f(E) =07(E) - o*(E); (8) 


then (7) can be rewritten 





A,(E - M/E) 
< 3)" insi(e- _)" eas 
_M? =f _E'ot(E’) | 
(E-E)s), 3 


"E 2 
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The condition (4) is clearly satisfied by the terms 


- containing D*(M) and the bound-state contribu- 


tions on the right-hand side of (9). This is also 
the case for the integral over o*: if in fact we 
subdivide the integration at an energy € suffi- 
ciently high so that for E’>«€, o*(E’) can be con- 
sidered a constant, then its contribution for 
E>e is essentially given by 


21 


red, > 


OME ME" == £o%(«), (10) 


Then, if we call 


cte)={" 
0 


condition (4) implies that 


lim EG(E) =constant. 
Ewa 


g(E')dE' 


(12) 


On the other hand, the following theorem on 
Stieltjes transforms holds*: 

Theorem: Let G(E) be the (convergent) integral 
defined by (11) (E,>0) and let 


lim zi, E'@(E")\dE' =0; (13) 
E-a 
then 
E 
lim (ce -f o(e as’) =0. (14) 
Eo E, 


The condition (13) is trivially satisfied in our 
case because of (2) [although it is less restric- 
tive than (2) since it implies only the existence 
of limits (1) and (2) in the mean]. From (14) 
and (12), the convergence of the integral 


(= -O(E') op, 
E 





eo 

f (EE (15) 
Ey 

follows, as was to be demonstrated.°® 

As a consequence, we see that the limiting 
behavior (6) is by no means compatible with the 
constancy of cross sections at high energies, 
ie., the difference o~(E) - ot(E) must go to zero 
faster than 1/InE. 

We want to note, finally, that an immediate 
consequence of our conclusion is that 


lim A(E) =0, 
E~o 


pb’ 


i) 


(16) 
where 
A(E) =[D(E) -D’(E)]/E. 


By using the unsubtracted dispersion relation 


(17) 





for D~(E) -D*(E) whose validity we have just 
discussed, we can write 
B, 2 —— p'f(E") 
AG) -2gt gts Jp BBP 


where B, are the constant residua of the bound- 
state contributions. Changing the integration 
variable to E”, the integral of the right-hand 
side of (18) is reduced to the form of a Hilbert 
transform. Its convergence and the possibility 
of inverting the Hilbert transform® are ensured 
by the existence of (15). This means that the 
expression 


(18) 


© A(E’) 
E” -E? 
must converge for all values of E (and go to 
zero for E~«), from which (16) follows. 

An empirical discussion on the limit value for 
A(E) was given by Goldberger et al.” in connec- 
tion with the possibility of using (18) as a sum 
rule. Here instead we have shown, from rather 
general theoretical ground, that (18) is valid 
without subtraction constants and that, as a con- 
sequence, the limit (16) holds. 


‘T. Ia. Pomeranchuk, J. Exptl. Theoret. Phys. 
(U.S.S.R.) 34, 725 (1958) [Soviet Phys. JETP 34(7), 
499 (1958)]. 

*Goldberger, Miyazawa, and Oehme, Phys. Rev. 99, 
986 (1955). See also: A. Salam,Proceedings of the 
CERN Symposium on High-Energy Accelerators and 
Pion Physics, Geneva, 1956 (European Organization 
of Nuclear Research, Geneva, 1956). Vol. II, p.179. 

3p, T. Matthews and A. Salam, Phys. Rev. 110, 
569 (1958); R. H. Dalitz, 1958 Annual International 
Conference on High-Energy Physics at Cern, edited 
by B. Ferretti (CERN Scientific Information Service, 
Geneva, 1958), p. 192; E. Galzenati and B. Vitale, 
Phys. Rev. 113, 1635 (1959); D. Amati, Phys. Rev. 
113, 1692 (1959). 

This is a slight modification of a thorem by G. H. 
Hardy and J. E. Littlewood, Proc. London Math. Soc. 
30, 34 (1930). Compare also with D. V. Widder, The 
Laplace Transform (Princeton University Press, 
1946), Lemma 5, p. 193. 

5(14) shows that if D*(E)/E satisfies only the weaker 
condition to be bounded at infinity, then the existence 
oi the integral (5) does not necessarily follow. The 
following example: o~ -o* ~ sin(InE)/InE (—0 for 
E-), for which EG(E) ~sin(InE), gives an oscillating 
behavior of f @(E')dE’ for E-«, On the other hand, 


if the limit (4) exists, then the integral (15) converges 
also in the more general case of bounded o* and o 
without the requirement of the actual existence of a 
limit for o~-0* atE—~«, This last fact is the con- 
sequence of another theorem of Hardy and Littlewood 
(reference 4). A trivial example: o~-o* ~sinE. 

SE, C,. Titchmarsh, Introduction to the theory of 
Fourier Integrals (Oxford Press, 1948), p. 724. 
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CAPTURE GAMMA RAYS FROM O** AND O** 
IN THE REGION OF THE GIANT RESONANCE. 
S. G. Cohen, P. S. Fisher, and E. K. Warburton 
[Phys. Rev. Letters 3, 433 (1959)]. 


On page 434, column 2, line 5, change 1:2 to 
2:1. 


EFFECTIVE EXCHANGE CONSTANT IN 
YTTRIUM IRON GARNET. D. T. Edmonds and 
R. G. Petersen [Phys. Rev. Letters 2, 499 
(1959). 


The value of 4nMg quoted below Eq. (6) should 
read 2449 gauss’ and not 1740 gauss. The value 
1740 gauss is the value of 47M¢ at room tem- 
perature rather than the desired value at 0°K. 
Making this correction, we have 


A(YIG) =0.270x10~* erg cm™, 
\ko! =5.46x10° cem7, 
A = 0.89 cm. 


We would like to thank Professor C. Kittel and 
Dr. W. P. Wolf who independently pointed out 
our error. 


'M. A. Gilleo and S, Geller, Phys. Rev. 110, 73 
(1958). 
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ABSTRACTS 








In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 

: published )”’ rather than to this Journal. 


IRREVERSIBLE THERMODYNAMICS AND CON- 

TINUUM MECHANICS. A. Cemal Eringen, Pur- 

due University, Lafayette, Indiana (Received 

June 11, 1959; revised manuscript received 
December 14, 1959). 


Principles of the irreversible thermodynamics 
; are recapitulated. Using the basic conservation 
principles and the principles of irreversible 
thermodynamics, the fundamental equations and 
constitutive relations are obtained for fluids, 
solids, and visco-elastic media. Equations for 
the entropy production and heat conduction are 
derived for various special media and for a 
general medium having internal constraints. For 
the latter case, a law generalizing Fourier’s 
law of heat conduction is found. 


OOP a 


ROLE OF MOLECULAR IONS, METASTABLE 
MOLECULES, AND RESONANCE RADIATION 
INTHE BREAKDOWN OF RARE GASES. A. V. 
Phelps, Westinghouse Research Laboratories, 
' Pittsburgh, Pennsylvania (Received August 14, 
1959). 





Theoretical calculations are made of the rates 
of arrival of atomic and molecular ions, of 
atomic and molecular metastables, and of non- 
resonance and resonance radiation at the cathode 
_ of parallel plane electrodes during the buildup 
' of current preceding the breakdown of a rare 

gas. The electron density is assumed to increase 
| exponentially with time and distance. Previous 
analyses are extended to obtain the particle cur- 
rents of molecular ions and of metastable atoms 
and molecules at the cathode. The Holstein- 
Biberman formulation of the transport equation 
for the density of atoms in the resonance state 
is solved numerically to find the fraction of 
resonance photons arriving at the cathode. The 
fraction of resonance radiation reaching the 
cathode varies with the growth constant of the 
electron current in very nearly the same manner 
as for delayed photons and varies much less 


OO 

















rapidly with electrode spacing and gas density 
than for the case of diffusing particles. The 
theory is applied to the calculation of the growth 
constant as a function of applied voltage for hel- 
ium at 100 mm Hg pressure and a 1-centimeter 
gap. The growth of current is found to be con- 
trolled by the rates of arrival of molecular 
metastables, molecular ions, and resonance 
radiation at the cathode. 


VAPOR PRESSURE OF He’ - He* SOLUTIONS. 
A. K. Sreedhar and J. G. Daunt, Department of 
Physics and Astronomy, Ohio State University, 
Columbus, Ohio (Received August 24, 1959). 


Measurements have been made of the vapor 
pressures of solutions of He*® in He* having molar 
He® concentrations 1.91, 3.34, 6.51, and 12.04% 
in the temperature range 1.4°K to 2.6°K. The 
measurements involved the observation of the 
excess pressure of the solution over that of pure 
liquid He* which was carefully maintained at the 
same temperature as the solution using an appa- 
ratus which ensured high precision of pressure 
transmission between the low- and high-tempera- 
ture parts of the apparatus, especially for liquid 
He* temperatures above the lambda temperature. 
The results are compared with previous evalua- 
tions and comparisons made with recent theories 
of He® - He* solutions. 


LONGITUDINAL HALL EFFECT. Ludwig 
Grabner, International Business Machines Re- 
search Laboratory, Poughkeepsie, New York 
(Received August 20, 1959). 


The “Hall field” is defined as the electric field 
which is an odd function of the magnetic field. 
Because of the Onsager relations it is perpen- 
dicular to the current. It is split into the con- 
ventional transverse Hall field (TH field) in the 
direction of B xI and a longitudinal Hall field 
(LH field) in the direction Ix(BxI). Some prop- 
erties of the LH field contrast with those of the 
TH field. In particular, it vanishes (a) at satura- 
tion, (b) for spherical energy surfaces, and (c) 
when the magnetic field is parallel to an axis of 
rotation of the crystal. Its dependence on car- 
rier concentration is the same as that of the TH 
field; its sign, however, is determined by band 
structure as well as carrier sign. For a system 
consisting of n-type germanium the properties 
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of the LH field are investigated theoretically 
and experimentally. 


F-CENTER GROWTH CURVES. P. V. Mitchell, 
D. A. Wiegand, and R. Smoluchowski, Carnegie 
Institute of Technology, Pittsburgh, Pennsylvania 
(Received June 17, 1959). 


It is shown that the increase in the density of 
F centers in KCl due to x-ray irradiation can 
be described kinetically by considering two in- 
dependent groups of negative-ion vacancies. One 
group is composed of vacancies in the lattice 
before irradiation, while the other is composed 
of vacancies generated during the irradiation. 
It is concluded from the difference in the elec- 
tron capture cross sections for the two groups 
of vacancies and the intensity dependence of the 
rate of generation of new vacancies that the 
vacancies are generated at dislocations by a 
multiple ionization process. 


SPECIFIC HEAT AND ELASTIC CONSTANTS 
OF CALCIUM FLUORIDE AT LOW TEMPERA- 
TURES. D. R. Huffman and M. H. Norwood, The 
Rice Institute, Houston, Texas (Received August 
20, 1959). 


Measurements of specific heat from 3.5°K and 
of elastic constants from 4°K to 300°K have been 
made on single crystals of calcium fluoride. 
The Debye characteristic temperature theta at 
0°K, 6,, as calculated from the specific heat 
data, is 508+5°K. Values of the elastic con- 
stants at 4.2°K are c,,=17.4, c,,=5.6, and c,, 
= 3.593 in units of 10°! dynes/cm*. 6, calculated 
from elastic constants is 513.6+2.5°K. The two 
6,’S are now in good agreement and are signifi- 
cantly greater than the value calculated from 
previous specific heat work. 


PROPERTIES OF CLEANED GERMANIUM SUR- 
FACES. R. Forman, Research Laboratories, 
National Carbon Company, Division of Union 
Carbide Corporation, Cleveland, Ohio (Received 
July 27, 1959). 


Field effect and conductivity measurements 
were made on cleaned (111) and (100) germanium 
surfaces. The surfaces were cleaned by an argon 
bombardment technique similar to that developed 


94 


—$—$$—__ 


by Farnsworth. The cleaned surfaces were ex. 
posed to a modified Brattain-Bardeen ambient — ra 
cycle and the change in surface properties 
measured. The cleaned surface data showed ter 
that a (100) surface is more highly p type than rol 
a (111) surface, and these results are interprety 7, 
qualitatively by an atomistic model. | wel 









FLUCTUATIONS IN CdS DUE TO SHALLOW | 
TRAPS. James J. Brophy and Robert J. Robin. | * 
son, Physics Division, Armour Research Mi 
Foundation, Chicago, Illinois (Received August | 

24, 1959). 


The observed current noise of a lightly doped RE 
cadmium sulfide single crystal under uniform RI 
illumination is quantitatively explained by appli- | No 
cation of the generation-recombination theorem | se 
to retrapping effects. The results yield a fre- Au 
quency factor of 10° sec™ for a discrete set of I 
traps located 0.4 ev below the conduction band ery 
and 10° sec™ for traps distributed between 0.3 30 
to 0.5 ev. When the quasi-Fermi level is near 


sul 
the discrete states, the current noise spectra col 
show a relaxation component characteristic of ma 
generation and recombination involving these oc 
levels. res 
ga] 


NUCLEAR MAGNETIC RESONANCE OF SOLID | 
HYDROGEN (67-86% ORTHO) AND SOLID DEv- | © 
TERIUM (33% AND 55% PARA). George W. duc 
Smith * and Robert M. Housley, The Rice Insti- 


tute, Houston, Texas (Received August 18, 960 
> AIC 
Nuclear magnetic resonance studies on solid BIL 


hydrogen (H,) of ortho concentrations from 67% | g, 
to 86% and of solid deuterium (D,) with para con- | py 
centrations of 33% and 55% are reported. The Cor 
method of Cunningham, Chapin, and Johnston Me 
was used to produce the high concentrations Aug 
ortho H, and para D,. The temperature (7) of 
appearance of side peaks associated with the T 
removal of rotational degeneracy in H, was found [| ¥4! 
to shift to 2.18°K for 86% ortho as compared to | ™ 
a value of 1.6°K for the normal concentration. fitt 
Our T, versus concentration data were found to Cu 
join smoothly to the previous lower concentration 
work of Sugawara et al. In addition an interesting f = 
hysteresis phenomenon was observed for H,. It 
was noted that the temperature of appearance of 
the side peaks when the H, is cooled is some 

0.15°K lower than the temperature of their dis- 
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appearance when the temperature is slowly 


’ raised. in the case of deuterium, no side peaks 


yere observed for either 33% or 55% para D, at 
temperatures down to 0.9°K. 7, was measured 
roughly for the two concentrations. For n-D,, 

T, was found to be 7.6+3 sec. This compares 
well with the theoretical value of 8 sec given by 
Moriya and Motizuki. For 55% para D,, a value 


of 3+1 sec was found. The theory of Moriya and 
Motizuki gives 5 sec. 
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*Now at General Motors Technical Center, Warren, 
Michigan. 


RECOMBINATION AND TRAPPING IN TELLU- 
RIUM. J. S. Blakemore, J. D. Heaps, K. C. 


_ Nomura, and L. P. Beardsley, Honeywell Re- 


search Center, Hopkins, Minnesota (Received 
August 24, 1959). 


Photoconductive decay in nearly perfect Te 
crystals shows a lifetime of up to 40 yusec at 
300°K. The temperature dependence does not 
support previous suggestions of radiative re- 
combination; at temperatures below 150°K 
marked trapping of excess minority electrons 
occurs, and at higher temperatures also the 
results suggest the activity of levels within the 
gap. Probably chemical impurities behave as 
recombination centers, and dislocations cer- 
tainly act in this capacity (lifetime is reduced 
to~1 psec when 10° dislocations/cm? are intro- 
duced). 


_ HIGH-TEMPERATURE MAGNETIC SUSCEPTI- 


BILITY OF SINTERED aMnS. John J. Banewicz, 
Southern Methodist University, Dallas, Texas, 
Robert Lindsay, Trinity College, Hartford, 
Connecticut, and Robert F. Heidelberg, Southern 
Methodist University, Dallas, Texas (Received 
August 20, 1959). 


The magnetic susceptibility of sintered aMnS 
was measured by a Faraday method at tempera- 
tures from 300°K to 1244°K. When the data were 
fitted by a least-squares method to the Weiss- 
Curie law, xq = Cyy/(T+ 6), using only points 
above 675°K, values of Cy, =3.94+0.02 and 
§=397+ 6°K were obtained. Even in this range 
the 1/y vs T curve was very slightly inflected 
towards the 1/y axis so the above values are 
considered only as establishing upper limits on 
the constants. It is suggested that the deviation 


of C M from the theoretical spin-only value of 
4.379 be interpreted in terms of a reduced ionic 
character of the MnS bond. 


WAVE FUNCTIONS FOR IRON d BAND. J. H. 
Wood, Solid-State and Molecular Theory Group, 
Massachusetts Institute of Technology, Cambridge, 
Massachusetts (Received August 28, 1959; re- 
vised manuscript received September 23, 1959). 


The one-electron wave functions for the d band 
of body-centered iron are examined and found to 
have significantly different character, depending 
on their energy location in the band. The bottom 
of the band is associated with diffuse wave func- 
tions, while the top is associated with more com- 
pact atomic-like functions. A smooth transition 
from one behavior to the other occurs as the 
energy is varied. This gives substantiation to 
recent descriptions of transition metals which 
have been based on “low d” assignments; here 
this is explained on the basis that a significant 
fraction of the d electrons in a solid do not “look” 
like atomic d electrons at all. 


THEORY OF MAGNETISM OF NiF,. Toru Mori- 
ya,* Bell Telephone Laboratories, Murray Hill, 
New Jersey (Received August 24, 1959). 


The magnetic properties of NiF, are studied 
theoretically with the use of a spin Hamiltonian 
approach. It is shown that the anisotropy of the 
susceptibility above the Néel temperature indi- 
cates a spin arrangement below the Néel tem- 
perature in which all the spins are perpendicular 
to the c axis and the angle between the two sub- 
lattice magnetizations is a little smaller than 7, 
giving rise to a small net ferromagnetic moment 
along the (100) direction. This model explains 
the torque data below the Neel temperature quan- 
titatively. The spin wave frequencies and the 
magnetic resonance frequencies are calculated. 
There is a low-frequency branch whose lowest 
frequency (k = 0) corresponds to the anisotropy 
energy in the ab plane. The magnetic resonance 
is expected at around 170 kMc/sec. A recent 
nuclear magnetic resonance measurement by 
Shulman seems to support the present model. 
The spin arrangement below the Néel tempera- 
ture proposed by Erickson from his neutron dif- 
fraction data seems to be neither possible theo- 
retically nor consistent with the other experi- 
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mental data. Possible structures of domains and 
domain walls are discussed, based on the same 
model. 


‘On leave of absence from Tokyo Metropolitan Uni- 
versity, Tokyo, Japan. 


OPTICAL CONSTANTS OF METALS. Maciej 
Suffczynski, * Department of Mathematics, 
Imperial College, London, England (Received 
August 17, 1959). 


A model is attempted for calculation of the 
interband contribution to the frequency-dependent 
dielectric constant of metals. The frequency 
region near the threshold for the interband 
transitions is considered. Emphasis in the 
model is laid on the bending of the energy bands 
near the Brillouin zone boundary. Attention is 
focused on cases when the Fermi surface ap- 
proaches the zone boundary or has a finite area 
of contact with it. The momentum matrix ele- 
ment is taken as constant, which is fitted so as 
to achieve agreement with the experimentally 
found dip in the dispersion curve of the extinc- 
tion coefficient. The values of the squared 
matrix element for the noble metals, copper, 
silver, and gold, which fit the experimental 
data of Schulz are found to be in the ratio 0.43: 
0.69:0.69. 


“Present address: Institute of Theoretical Physics, 
University of Warsaw, Warsaw, Poland. 


SUPERCONDUCTIVITY OF CONTACTS WITH 
INTERPOSED BARRIERS. Hans Meissner,* 
Department of Physics, The Johns Hopkins Uni- 
versity, Baltimore, Maryland (Received 
August 25, 1959). 


Resistance vs current diagrams and “diagrams 
of state” have been obtained for 63 contacts be- 
tween crossed wires of tin. The wires were 
plated with various thicknesses of the following 
metals: copper, silver, gold, chromium, iron, 
cobalt, nickel, and platinum. The contacts be- 
came superconducting, or showed a noticeable 
decrease of their resistance at lower tempera- 
tures if the plated films were not too thick. The 
limiting thicknesses were about 35 x10~* cm for 
Cu, Ag, and Au; 7.5x10~* cm for Pt, 4x10 
cm for Cr, and less than 2x10~° cm for the 
ferromagnetic metals Fe, Co, and Ni. The in- 
vestigation was extended to measurements of 
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the resistance of contacts between crossed wires 
of copper or gold plated with various thicknesses 
of tin. Simultaneous measurements of the (lon- 
gitudinal) resistance of the tin-plated gold or 
copper wires showed that these thin films of tin 
do not become superconducting for thicknesses 
below certain minimum values. These latter 
findings are in agreement with previous meas- 
urements at Toronto. The measurements at 
Toronto usually were believed to be unreliable 
because films of tin evaporated onto quartz sub- 
strates can be superconducting at thicknesses 
as small as 1.6x10~° cm. It is now believed 
that just as superconducting electrons can drift 
into an adjoining normal conducting layer and 
make it superconducting, normal electrons can 
drift into an adjoining superconducting layer 
and prevent superconductivity. 


‘Now at Department of Physics, Stevens Institute of 
Technology, Hoboken, New Jersey. 


EXCHANGE EFFECTS IN ALKALI OZONATES. 
S. A. Marshall and D. O. Van Ostenburg,* Phys- 
ics Division, Armour Research Foundation, 
Chicago, Illinois (Received August 24, 1959). 


Resonance absorption has been observed in 
polycrystalline samples of potassium and cesium 
ozonate. The absorptions are characterized by 
nearly free electron spin g values with line widths 
ranging from 26 to 49 oersteds. Analysis of line 
shapes based on the Anderson-Weiss model for 
exchange interaction between equivalent mag- 
netic ions and calculations of line-shape moment 
ratios based on the Van Vleck theory, strongly 
suggest exchange-interaction narrowing. Values 
of exchange frequencies for potassium and 
cesium ozonate are estimated to be 6x10" rad 
sec™ and 3x10" rad sec”, respectively. 

“Present address: Argonne National Laboratory, 
Lemont, Illinois. 


ASPECTS OF THE THEORIES OF DISLOCATION 
MOBILITY AND INTERNAL FRICTION. Jens 
Lothe, Institute of Physics, University of Oslo, 
Oslo, Norway (Received August 28, 1959). 


The theory of dislocation mobility is recon- 
sidered, and it is concluded that the interaction 
between thermal waves and a moving vibrating 
dislocation causes a drag to the first order in 
v/c. Modifications of the Seeger-Donth theory 
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for the Bordoni peak are suggested. When 
account is taken of the diffusion of kinks, general 
agreement with experiment is obtained. The 
theory of internal friction in the megacycle 

region is briefly reviewed and discussed. No 
thorough comparison with experimental data 

has been performed. 


DYNAMICAL THEORY OF DIFFUSION IN CRYS- 
TALS. Oscar P. Manley, Energy Conversion 
Group, Massachusetts Institute of Technology, 
Cambridge, Massachusetts, and Stuart A. Rice, 
Department of Chemistry and Institute for the 
Study of Metals, University of Chicago, Chicago, 
llinois (Received August 27, 1959). 


Corrections to a previous paper by the first 
author are presented. A calculation leading to 
the simplification of the pair correlation functions 
is also presented. 


THEORY OF FERROMAGNETIC RESONANCE 
INRARE EARTH GARNETS. III. GIANT ANI- 
SOTROPY ANOMALIES. C. Kittel, Miller Insti- 
tute for Basic Research in Science, University 
of California, Berkeley, California (Received 
August 17, 1959). 


The giant sharp anisotropy anomalies in yttrium 
iron garnet discovered by Dillon and attributed to 
rare earth ions are explained in terms of the cross- 
oers or, most probably, near crossovers of the 
energy levels of the rare earth ions in the com- 
bined crystal and exchange fields. The central 
consequence deduced by the theory, independent 
of the detailed behavior of the energy levels, is 
that -H (A@y2)* =2£kTN’/M, where H , is the 
anisotropy field; Aé,,. the angular width of the 
peak; N’ the concentration of rare earth ions at 
the crossover; £=1 for an isothermal process 
and $ for an adiabatic process; M is the total 
magnetization. The result assumes kT > level 
splitting at apparent crossover. This and other 
results appear to be in satisfactory agreement 
vith the available experimental data. Crystals 
vith giant anomalies have useful properties for 
uiabatic demagnetization experiments —the 
ferric-rare earth exchange field of 100 koe can 
ve effectively turned on or off by rotation of the 
magnetic moment of the crystal. In an Appendix 
We examine the validity of the molecular field 
approximation to the ferric-rare earth ion inter- 








action, and we find that the approximation is ex- 
cellent. A further Appendix discusses the ques- 
tion of crossing or noncrossing of energy levels 
in static fields. 


ELECTRON-SPIN PARAMAGNETIC RESONANCE 
STUDIES OF METALLIC LITHIUM IN NEUTRON- 
IRRADIATED LiF. Y. W. Kim, R. Kaplan, and 
P. J. Bray, Department of Physics, Brown Uni- 
versity, Providence, Rhode Island (Received 
August 31, 1959). 


The electron-spin paramagnetic resonance ab- 
sorption of conduction electrons in lithium metal 
in neutron-irradiated LiF has been observed. 
The dependence of the electron-spin paramag- 
netic resonance line shape and width on the tem- 
perature during irradiation and subsequent ther- 
mal annealing is discussed. An upper bound is 
set on the probability for flip of the spin in a 
collision with the metal surface. 


FERROMAGNETIC ANISOTROPY IN CUBIC 
CRYSTALS. F. Keffer and T. Oguchi,* Sarah 
Mellon Scaife Radiation Laboratory, University 
of Pittsburgh, Pittsburgh, Pennsylvania (Re- 
ceived August 21, 1959). 


A re-evaluation has been made of Van Vleck’s 
second-order perturbation theory of dipolar- 
type anisotropy in cubic ferromagnets. In the 
low-temperature limit of strong correlation be- 
tween the direction of neighbor spins, the first 
anisotropy constant K, varies as the tenth power 
of the magnetization. The theory is somewhat 
analogous to a previous treatment of quadrupolar- 
type anisotropy in the strong-correlation limit. 
In both cases, the results are in agreement with 
the Akulov-Zener classical theory. For the 
dipolar case, complete agreement is also esta- 
blished between the Dyson-type spin-wave anal- 
ysis of Charap and Weiss and the Holstein- 
Primakoff approach. Higher order terms in the 
latter are shown to lead to the Charap-Weiss 
correction from exchange interaction between 
spin waves, and this correction is extended to 
S>1/2. Essentially the same correction is ob- 
tained very easily from a simple modification of 
the Van Vleck formalism to take careful account 
of the average energy involved in simultaneous 
reversal of neighbor spins. It is shown that 
spin-wave theory, in agreement with classical 
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theory, predicts identical values of dipolar-type 
anisotropy whether measured statically in a 
torque experiment or dynamically in a micro- 
wave resonance experiment. 


"Permanent address: Department of Physics, Tokyo 
University of Education, Ohtsuka-kubo-machi, Bunkyo- 
ku, Tokyo, Japan. 


MAGNETIC - FIELD DEPENDENCE OF THE 
ULTRASONIC ATTENUATION IN METALS. 
Morrel H. Cohen and Michael J. Harrison, Insti- 
tute for the Study of Metals, University of Chi- 
cago, Chicago, Illinois, and Walter A. Harrison, 
General Electric Research Laboratory, Schenec- 
tady, New York (Received August 31, 1959). 


A self-consistent, semiclassical treatment is 
given for the attenuation of a sound wave by a 
free-electron gas in a positive background which 
supports the sound wave. Emphasis is placed 
upon the kinds of magnetic-field dependence 
which can be found under a wide range of values 
for the magnetic field, frequency, and mean free 
path. Applications of the general formalism 
include propagation parallel and perpendicular 
to the magnetic field. The special phenomena 
studied include geometric resonances, cyclotron 
resonances, and magneto-plasma resonance. A 
qualitative physical interpretation of the various 
effects found in the detailed calculations is also 
presented. 


PHOTOCONDUCTIVITY OF CUPROUS OXIDE IN 
RELATION TO ITS OTHER SEMICONDUCTING 
PROPERTIES. Frank L. Weichman,* North- 
western University, Evanston, Illinois (Received 
May 20, 1959). 


It is shown that at least two of the three known 
peaks in the photoconductivity versus wavelength 
curve of Cu,O are strongly dependent on the oxy- 
gen content of the material. The peak near 600 
my can be shifted from shorter to longer wave- 
lengths by exposing the sample to relatively low 
oxygen pressures at 860°C. The peak near 800 
muy is found to increase in magnitude relative to 
the others under these circumstances. This is 
in apparent contrast to the results of previous 
investigators. It is further shown that both re- 
sults are consistent with the values of the high- 
temperature activation energy of the conductivity 
process, as well as with the change in bandgap 
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as a function of temperature as obtained from 
optical absorption measurements under the as- 
sumption that Cu,O can be classified as a nonpoly 
semiconductor. 


‘Now at the Department of Physics, University of 
Alberta, Edmonton, Alberta, Canada. 


HARTREE-FOCK CALCULATIONS FOR Mn*+ 
IN CUBIC FIELDS. R. E. Watson, Department o 
Physics, Massachusetts Institute of Technology, 
Cambridge, Massachusetts (Received August 31, 
1959). 


Hartree-Fock calculations have been carried 
out for the Mn** ion in cubic fields produced 
by sets of octahedrally placed point charges. 
The two types of cubic 3d electrons were allowed 
to have different radial dependence. The relax- 
ation of the constraint of having a single 3d radial 
function led to effects which were small but which 
the writer believes should be included in any de- 
tailed treatment of a paramagnetic ion in a crys- 
talline field. 


CHARGE TRANSFER IN MOLECULAR HYDRO- 
GEN. T. F. Tuan,” University of Pittsburgh, 
Pittsburgh, Pennsylvania, and E. Gerjuoy,! 
General Atomic, San Diego, California (Re- 
ceived August 17, 1959). 


The effects of the molecule on electron capture 
by protons in hydrogen gas have been investigated 
in first Born approximation with different types 
of electronic wave functions. It always has been 
supposed that if the incident proton velocity is 
large compared to electronic velocities molecular 
effects may be neglected, and that one may then 
assume one H, molecule is equivalent to two 
hydrogen atoms for purposes of charge transfer. 
Instead it appears that charge transfer in H, at 
high energies bears no simple relationship to 
charge transfer in atomic hydrogen. In particular, 
among other effects: (i) in the high-energy limit 
20 =1.2 - 1.40 4; (ii) at lower energies there is 
important interference between the capture am- 
plitudes from the two atoms in the molecule. It 
is also found that transitions to ungerade states 
of H,*, although unimportant in the energy range 
of present experiments, become appreciable at 
high energies. 


* 
Present address: Physics Department, Northwestern 
University, Evanston, Illinois. 
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ton leave from the University of Pittsburgh, Pitts- 
burgh, Pennsylvania. 


DSTORTED-WAVE METHOD FOR ELECTRON 
CAPTURE FROM ATOMIC HYDROGEN. R. H. 
Bassel, University of Pittsburgh, Pittsburgh, 
Pennsylvania, and E. Gerjuoy, f General Atomic, 


The distorted-wave method, in which distorted 
waves replace plane waves in the matrix elements 
for the collision amplitude, is given a rigorous 
formulation for arbitrary rearrangement colli- 
sions. In atomic collisions involving electron 
rearrangement, the formulation enables removal 
of the unphysical internuclear Coulomb interaction 
from the “perturbation” Hamiltonian. The pro- 
cedure is applied to the problem of electron cap- 
ture by protons in atomic hydrogen, eliminating 
explicit dependence on the proton-proton inter - 
action. Numerical results of a first order cal- 
culation are in good agreement with experiments 
from 35 to 200 kev. Despite this agreement and 
the agreement of previous calculations it is felt 
that the charge transfer process is still not well. 
understood. 

‘Present address: Electronuclear Research Division, 
Oak Ridge National Laboratory, Oak Ridge, Tennessee. 
ton leave from the University of Pittsburgh, Pitts- 

burgh, Pennsylvania. 


ATOMIC PHOTOELECTRIC EFFECT AT HIGH 
ENERGIES. R. H. Pratt,* The Enrico Fermi 
Institute for Nuclear Studies and Department of 
Physics, The University of Chicago, Chicago, 
llinois (Received June 29, 1959; revised manu- 
script received December 28, 1959). 


Total cross sections are obtained for the photo- 
electric effect from the K shell of an atom of 
arbitrary charge, in the limit of high energies. 

An approximate analytic formula then is deduced 
to cover the entire high-energy region. For 

heavy elements and very high energies the dif- 
ferences from previous predictions are large. 

tis noted that these results also apply to other 
Processes, including the one-photon annihilation 
of fast positrons. 








tern 








* 
Now at Department of Physics, Stanford University, 
Stanford, California. 


San Diego, California (Received August 17, 1959). 





ISOTOPE SHIFT IN THE Mo SPECTRUM. Hack 
Arroe and John M. Cornwall, Denver Research 

Institute, University of Denver, Denver, Colo- 

rado (Received August 26, 1959). 


The isotope shift in the molybdenum spectrum 
was measured with a Fabry-Perot interfero- 
meter, using separated isotopes. For the line 
16678 cm™~', a definite minimum isotope shift of 
-8.9x107* cm~! was observed between Mo™ and 
Mo”. The negative sign indicates that the direc- 
tion of the shift is opposite to that predicted by 
the volume effect. 


CONTRIBUTION TO THE ELECTROSTATIC 
SELF-ENERGY OF A CHARGED LIQUID DROP. 
W. D. Foland, University of Massachusetts, 
Amherst, Massachusetts (Received August 27, 
1959). 


The electrostatic self-energy of the liquid- 
drop-model nucleus has been divided into two 
parts, one of which has been evaluated. The 
problems arising due to difficulties in choosing 
limits of integration—difficulties which occur due 
to the presence in the integrand of the factor 
Y,,~' —are all contained in the formulation of 
that part of the self-energy which we have ob- 
tained. The evaluation of the other part of the 
self-energy poses no such limit problems. 

The shapes which we consider for the drop 
are the shapes adopted by a nucleus undergoing 
symmetric fission. 

The self-energy is expressed as a multiple- 
power series in the deformation parameters a,. 
The part of the self-energy which we calculate 
is zero through the general seventh-power term 
in the multiple-power series for the self-energy. 


ANOMALOUS ELECTRIC DIPOLE CONVERSION 
COEFFICIENTS IN ODD-MASS ISOTOPES OF 
THE HEAVY ELEMENTS. F. Asaro, F. S. 
Stephens, J. M. Hollander, and I. Perlman, De- 
partment of Chemistry and Lawrence Radiation 
Laboratory, University of California, Berkeley, 
California (Received August 13, 1959). 


A detailed review is given of experimental 
data on the anomalous L- and M-shell conver - 
sion coefficients of low-energy electric dipole 
transitions observed in the decays of odd-A 
nuclei of high atomic number. 
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The data are consistent in every case with the 
interpretation that the £1 conversion coefficients 
in the Ly shell agree with the theoretical, 
model-independent coefficients by Sliv and Band 
and by Rose. It is definitely established in 
several well-measured cases that the Ly and Lyy 
coefficients are substantially larger than the 
theoretical values. The most striking anomaly 
occurs in the 84.2-kev transition in Pa’, where 
the Ly and Ly coefficients are 21 and 15 times 
larger than the theoretical values, respectively. 
The experimental L; and Ly coefficients are 
correlated with the lifetimes of the transitions, 
and it is shown that the magnitude of the anomaly 
(Ly plus Ly) is proportional to the retardation 
in gamma-ray lifetime over that calculated from 
the single-proton formulas. No systematic 
trend has been observed in the deviations of the 
Ly and Ly coefficients individually. 


STUDY OF SOME REACTIONS OF 14.8-Mev 
DEUTERONS WITH THE LITHIUM ISOTOPES. 
E. W. Hamburger and J. R. Cameron,* Univer- 
sity of Pittsburgh, Pittsburgh, Pennsylvania 
(Received August 25, 1959). 


Metallic targets of enriched Li® and Li’ were 
bombarded with 14.8-Mev deuterons. The out- 
going particles were magnetically analyzed and 
detected in nuclear emulsions or scintillation 
counters. Angular distributions for the following 
reactions were measured: Li®(d,p)Li™*, 4.6 and 
7.46 Mev; Li’(d,p)Li®*, 2.28 Mev; Li"(d, t)Li™, 
3.57 Mev; Li®(d,d’)Li®™*, 2.19 Mev; Li®(d, a)He*, 
ground state; Li®(d, ¢)Li®, and Li®(d,He*)He®. A 
search was made for other levels produced in 
these reactions and in the Li"(d,d’)Li™* reaction. 
No new levels were observed; upper limits on 
the cross sections of unobserved levels are 
given. The results of the (d,p), (d,t), and (d, He*) 
reactions were compared with stripping theory 
and with the predictions of the shell model. The 
agreement is, in general, satisfactory. 

The Li®(d, ¢) and Li®(d, He*) reactions produce 
a continuous spectrum of outgoing particles at 
each angle. The angular distribution of the tri- 
tons (or He® particles) is strongly peaked for- 
ward and can be fitted with Butler curves having 
l=1, r5=6x10-** cm. The spectra at small 
angles show a peak corresponding to the Li® (or 
He*) ground state, but at large angles the peak 
disappears and the spectra are constant. The 
shape of the spectra can be explained by the 
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a 


Butler formula if the energy dependence of the 
reduced width is taken into account. The form 
of this dependence was deduced from the results !€¥ 
of the nucleon-helium elastic scattering experi- Er’ 
ments. Er 





"Now at the University of Wisconsin, Madison, Wig. (ga 
consin. 


RADIATION FROM Ga™. Thomas H. Jacobi, * 
Department of Physics, University of California, 
Los Angeles, California, H. A. Howe, U. S. PH 
Naval Radiological Defense Laboratory, San : 
Francisco, California, and J. Reginald Richar 
Department of Physics, University of California, ' Nu 
Los Angeles, California (Received July 6, 1959), Ch 


The 8* decay of Ga® and the subsequent gamm 
spectrum from the de-excitation of Zn™ have 
been studied. The 8* transition to the ground 
state of Zn™ has an end-point energy of 6.05+0.03. ph 
Mev with log ft =6.6 and another 8* group or Me 
groups has an end-point energy of 2.79+0.08 Mey sp 
with a log ft=4.6. No 8* transition to the first | yj 
excited state of Zn(2+) was observed and it was ys' 
shown that the log ft corresponding to this tran- | a¢ 
sition is greater than 7.7. The gamma rays that | to 
were observed (0.80, 0.99, 1.25, 1.38, 1.56, us 
1.78, 1.95, 2.18, 2.34, 2.99, and 3.32 Mev) in- 
dicate excited states at 0.99(2+), 1.78(2+), and — th 


yr 


oS 
bi] 





2.35(2+) Mev. gu 
The positron results indicate that the ground yj 
state of Ga“ has spin and parity 0+. an 


_ 
Now at Emerson-Electric of Saint Louis, Saint Louis, 
Missouri. 


ELECTRON CAPTURE DECAY OF Tm’® AND - 
Tm", K. P. Jacob and J. W. Mihelich, Univer- ™ 
sity of Notre Dame, Notre Dame, Indiana, and. | 
B. Harmatz and T. H. Handley, Oak Ridge Na- m 
tional Laboratory, Oak Ridge, Tennessee (Re- _{o 
ceived August 19, 1959). 1 


The electron capture of Tm’® (87-day) to levels | ™ 
in Er’® was investigated with permanent magnet « 
spectrographs and coincidence scintillation spec- 
trometers. The following levels in Er’® have 
been established: 79.8 (2+), 264.3 (4+), 548.9 
(6+), 822.4 (2+), 897.0 (3+), 996.2 (4+), 1095.1 
(3- and T,,.=1.2x10~ sec) and 1543.1 (3-) kev. 
The internal conversion data for Tm!® (7.7-hour) 
suggest levels in Er’ at 80.6 (2+), 265.1 (4+), 
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545.3 (6+), 787.1 (2+), 860.6 (3+), and 957.2 (4+) 
kev, with many more high-lying levels. Energy 
level schemes are proposed for both Er’® and 
gr, The levels at 822, 897, and 996 kev in 
er’ and 787, 861, and 957 kev in Er’® may 
possibly be associated with electric quadrupole 
 Wis- (gamma) vibrations. Some general features re- 
garding these vibrational levels are discussed 
and compared with available data on other even- 
‘ even nuclei in the rare earth region. 


= 
the 
orm 
eSults 
peri- 


S. | pHOTODISINTEGRATION OF THE DEUTERON 
nN | wirH 94- Mev BREMSSTRAHLUNG RADIATION. 
hards y 4. Galey, The Enrico Fermi Institute for 
ornia, " yyclear Studies, The University of Chicago, 
1959), Chicago, linois (Received August 24, 1959). 


amma Differential cross sections for the reaction 

re y+d~p+n have been determined at laboratory 

nd angles of 45°, 75°, 90°, and 135° for laboratory 

D + 0.03 photon energies from approximately 50 to 90 

P Mev. At each of the above angles the energy 

8 Mev | spectrum of the recoil protons was determined 

rst | with a counter telescope and a pulse-height anal- 

was ysis system. The low cross section for this re- 

ran- | action necessitates the reduction of background 

that § toa minimum. This was accomplished by the 

y use of a gaseous target and a particle selection 

in- ‘technique. The differential cross sections and 

ind _ the estimates of the parameters describing an- 
‘gular distributions are in reasonable agreement 

ind with recent calculations by de Swart and Marshak 
and by Zernik, Rustgi, and Breit. 


| Louis, 





NUCLEAR STRUCTURE EFFECTS IN T!?°, 

B. I. Deutch* and N. Goldberg, Physics Depart- 
ID ment, University of Pennsylvania, Philadelphia, 
ver- | Pennsylvania (Received August 24, 1959). 


nd. Angular correlation measurements have been 
a- made on the 404-kev and 279-kev transitions 
€- following electron capture from 52-hour Pb?®. 


| The following directional correlations were meas- 
evels | ed: gamma-gamma, K conversion electron - 
snet « gamma, and gamma -A& conversion electron. The 
pec- angular distributions of all measurements were of 

the form W(@)=1+AP,(cos@). The following 

| results were obtained. For the 404, -279Ke™ 
1 correlation A= -0.052+0.015, for the 279y - 404Ke- 
correlation A= -0.036+ 0.010, and for the 404, 
-279 correlation A= -0.151+0.010. No effects 
tue to perturbations of the intermediate state 


V. 
our) 


), 










were found. 

Using the well- measured 279-kev K conversion 
coefficient, the 404-kev K conversion coefficient 
was determined to be 0.1174 0.015 while the 680- 
kev K conversion coefficient was determined to 
be 0.011+0.004. From the gamma-gamma 
angular correlation measurement results and 
the mixture ratio in the 279-kev transition deter - 
mined by Stelson and McGowan, the mixture 
ratio in the 404-kev transition was determined 
to be 6, =+0.043+0.010. The Sliv value for the 
K conversion coefficient for this mixture ratio 
is 0.147. The correlations involving conversion 
electrons rule out a 3/2 - 3/2 ~—1/2 transition. 
They indicate that the M1 particle parameter in 
the 404-kev transition is about twice the theoret- 
ical prediction. The particle parameters for 
the 279-kev transition also do not agree with 
the theoretical predictions. 


“Present address: Bartol Research Foundation, 
Swarthmore, Pennsylvania. 


RADIOACTIVE DECAY OF LUTETIUM 173. 

R. G. Wilson and M. L. Pool, Department of 
Physics and Astronomy, Ohio State University, 
Columbus, Ohio (Received August 24, 1959; re- 
vised manuscript received October 12, 1959). 


Ytterbium oxide enriched to 92.6% in the 173 
mass number was irradiated with 6-Mev protons. 
An activity with a half-life of 625+ 50 days was 
produced and its assignment to Lu’”® confirmed 
by the identification of the ytterbium K x-ray, 
gamma rays corresponding to transitions among 
the energy levels established in Yb*** by Coulomb 
excitation of the separated isotope, and the activ- 
ities produced by the similar proton irradiations 
of the other enriched isotopes of ytterbium. The 
observed activity of Lu’”* consists of the L and 
K x-rays of ytterbium and gamma rays with en- 
ergies of 7941, 10141, 17242, 273+2, 34924, 
~450, 556+4, and 630+ 4 kev and two other 
gamma rays not resolved in the gamma-ray 
spectrum but observed in gamma-gamma coin- 
cidence measurements. These two gamma rays 
have energies of approximately 179 and 282 kev. 
Because no positron radiation exists in the activ- 
ity of Lu’”*, the mode of decay is solely by elec- 
tron capture to Yb’. Gamma-gamma coinci- 
dence measurements have lead to the assignment 
of a 633-kev level with a probable spin of 5/2+ 
and the confirmation of a 351-kev level in Yb'”* 
in addition to the previously known 179.5- and 
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78.7-kev levels. The disintegrations of Lu’” to 
the 633-kev level are by L capture only. Branch- 
ing ratios for the electron capture transitions of 
Lu?” to the levels of Yb’”* and approximate rela- 
tive probabilities for the transitions in Yb*”* are 
given in a proposed energy level scheme for the 
decay of Lu’”*. The choice of 9/2- for the 
ground-state spin of Lu’” is the most consistent 
with the proposed energy level scheme. 


CONTINUOUS ELECTRON SPECTRUM ACCOM- 
PANYING K CAPTURE. H. Daniel,* G. Schupp, 
and E. N. Jensen, Institute for Atomic Research 
and Department of Physics, Iowa State Univer- 
sity of Science and Technology, Ames, Iowa 
(Received August 31, 1959). 


The continuous electron spectrum accompany- 
ing K capture has been observed, in the case of 
Cs'', in three different ways: by measuring 
coincidences between one K x ray and an elec- 
tron, by measuring coincidences between two K 
xX rays and an electron, and by measuring the 
single-electron spectrum. Furthermore, coinci- 
dences between two K x rays have been meas- 
ured. The coincidence data were taken with 
scintillation counters while the single-electron 
spectrum was taken with a magnetic spectrom- 
eter. The electron spectrum, including the ab- 
solute intensity, and the probability for a double 
hole in the K shell are both found to be in agree- 
ment with the theory of Primakoff and Porter. 


* 
On leave from Max Planck Institute for Nuclear 
Physics, Heidelberg, Germany. 


SCATTERING OF 43-Mev a PARTICLES BY 
NUCLEI J. L. Yntema, B. Zeidman, and B. J. 
Raz,* Argonne National Laboratory, Lemont, 
Illinois (Received August 5, 1959). 


The differential scattering cross sections for 
the elastic scattering of 43-Mev a particles by 
Zn, Ag, Rh, and Au are obtained over the angu- 
lar range from 12° to 55°. The absolute errors 
in the cross sections are estimated to be less 
than 2% over most of the angular range. The 
group structure in the inelastic spectrum ob- 
served with protons and deuterons is also ob- 
served with a particles. In the spectra of Ag, 
Rh, and Zn one group near 2.7 Mev was observed. 
In Al, V, and Ni a number of groups were ob- 
served. In Au and Ta no groups could be re- 
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solved. Angular distributions were obtained for 
the 2.5-Mev group in Rh and Ag and the 2.8-Mey 
group in Zn. In the case of the 2.8-Mev group ip 
Zn, the angular distribution shows a very marke} 
diffraction pattern. The data were fitted with 
Butler’s theoretical expression modified to elim. 
inate the assumption that the scattering of 43- 
Mev a particles by nucleons is isotropic in the 
center-of-mass system. It is shown that this 
assumption can be easily eliminated in the for- 
mal treatment. A comparison with the angular 
distributions obtained for the 1.04-Mev 2* leve| 
in Zn with both 43-Mev a particles and 21.6-Mey 
deuterons indicates that /=3 is the most plausib|: 
assignment for the 2.8-Mev group. 


* ' 
Permanent address: State University, College on 
Long Island, Oyster Bay, New York. 


p’-y ANGULAR CORRELATION IN THE INELAS. 
TIC SCATTERING OF 16.6-Mev PROTONS BY 
Mg”*. Hajimé Yoshiki,* Palmer Physical Labo- 
ratory, Princeton University, Princeton, New 
Jersey (Received August 10, 1959). 


The angular correlation between the inelas- 
tically scattered protons leaving the first excited 
state of Mg™* (1.37 Mev, 2+) and its de-excitation 
gamma rays produced by the bombardment of 
16.6- Mev protons on a natural magnesium target 
has been investigated. For 6»’=30°, 42.5°, and 
70° the result shows that the correlation functions 
are represented fairly well by the simple Born 
approximation calculation, sin?[2(6, - ¢,)|. 
However, at larger angles, 6’=95°, 120°, and 
150°, changes in the shapes of the correlation 
curves are observed, introducing a contribution 
from a sin*(6, -6,') term. This fact suggests 
the existence of exchange effects as well as the 
distortion effect in the present nuclear reactions. 

"Present address: Pupin Physics Laboratory, Colum-, 
bia University, New York, New York. 


NUCLEAR REACTIONS OF HIGH-ENERGY C”, 
N**, AND O'* WITH CARBON. Murry A. Tamers, 
Department of Chemistry, University of Texas, « 
Austin, Texas, and Richard Wolfgang, Depart- 
ment of Chemistry, Yale University, New Haven, 
Connecticut (Received August 25, 1959). 


The reactions of 160-Mev O'*, 140-Mev N*, 
and 120-Mev C™ with C” to yield radioactive 
products have been studied and the results com- 
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pared with those for proton- induced reactions 
giving formally similar compound systems. Ex- 
citation functions for products formed at lower 
energies are similar in the corresponding pro- 
ton and heavy-ion cases, indicating that their 
formation involved similar low-energy, low-spin 
compound nuclei. However, no “tail” due to 
knockon cascade processes was observed with 
heavy ions. Yields for higher energy products, 


- chiefly F**, were much higher in the heavy-ion 


bombardments. This may be attributed partly to 
higher average energy deposition with heavy 
ions, due to the absence of knockon processes, 
and partly to enhanced alpha emission from the 
distorted high-spin compound nuclei formed by 


, heavy ions. The contribution of various stripping 
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processes to these yields is discussed. 


EXCITATION FUNCTIONS FROM THE ALPHA- 
PARTICLE BOMBARDMENT OF TITANIUM -46. 
D. O. Raleigh* and J. M. Miller, Chemistry 
Department, Columbia University, New York, 
New York (Received August 17, 1959). 


Excitation functions from a-particle bombard- 
ment of Ti*® were obtained for the (a, p), (a,n), 
(a,pn), (@, 2m), (a@,p2n), and (a, apn) reactions 
from 10 to 40 Mev a-energy, using the stacked- 
foil method. The results indicate that evapora- 
tion theory supplies an accurate picture of neu- 
tron-proton competition in nuclear evaporation, 
but that direct-interaction processes can com- 
pete significantly with compound-nucleus forma- 
tion. 


*Now at Atomics International, Canoga Park, Cali- 
fornia. 


NUCLEAR SPIN RELAXATION AND NUCLEAR 
ELECTRIC DIPOLE MOMENTS. E. M. Purcell, 
Harvard University, Cambridge, Massachusetts 
(Received August 26, 1959). 


Proposals that nuclear spin relaxation in an 
appropriate system could serve as a test for the 
existence of a nuclear electric dipole moment 
are examined with attention to the consequences 
of the fact that the electric field at the nucleus is 
proportional to the nuclear acceleration. It is 
found that low-frequency fluctuations of the local 
electric field are suppressed. In particular, the 


necessarily negative correlation of the momentum 


transferred in consecutive collisions of an atom 


in a gas alters the spectral density of the pertur- 
bation, from that of uncorrelated pulses, by the 
factor w*7,’/(1+w*7,”), where 7, is the mean 
time between collisions. It follows that fairly 

low gas density is preferable to high. At optimum 
density a light gas at room temperature carrying ' 
nuclear electric dipole moments of magnitude 

e x10-** cm should have a spin relaxation time, 

in the absence of competing processes, around 

10 minutes. A formula is given for the electri- 
cally induced spin relaxation rate in a crystal. 
The process is hopelessly slow. In the electric 
coupling of the lattice vibrations to the spin, the 
ordinarily dominant “two-phonon” or “Raman” 
process is absent, because of the linearity of the 
connection between local electric field and nu- 
clear motion. 


ISOMERIC STATES OF Nd"! AND Sm’. RalphA. 
James* and Carleton D. Bingham, Department of 
Chemistry, University of California, Los Ange- 
les, California (Received August 17, 1959). 


Two previously unreported nuclear isomers 
which belong to the series having 81 neutrons 
and even numbers of protons have been found and 
studied. These are Nd’*” (half-life 63.7 sec- 
onds, gamma-ray energy 0.76 Mev) and Sm**5”” 
(half-life 137 seconds, gamma-ray energy 0.68 
Mev). The observed properties agree with the 
assignment of h,,,. and dy, to the upper and lower 
states, respectively, as predicted by the single- 
particle model. 


* 
Permanent address: Lawrence Radiation Laboratory, 
Livermore, California 


EFFECT OF COMPRESSION ON THE DECAY 
RATE OF Tc” METAL. R. A. Porter, Theo- 
retical Division, Lawrence Radiation Labora- 
tory, University of California, Livermore, Cal- 
ifornia, and W. G. McMillan, Department of 
Chemistry, University of California, Los Ange- 
les,and the RAND Corporation, Santa Monica, 
California (Received August 28, 1959). 


Motivated by Bainbridge’s measurement, we 
have attempted a theoretical calculation of the 
change in lifetime of the internal conversion of 
Tc*™” in compressed Tc metal. We have em- 
ployed the Thomas- Fermi statistical potential, 
corrected for the self-potential of the electron 
in question, to obtain the initial and final state 
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electronic wave functions for two volumes: the 
normal uncompressed state, and 10% compres- 
sion. Because the energy available is so low 
(~2 kev), the only contributions to the internal 
conversion coefficient come from the M and 
higher shells, and mainly from the 3p and 3d 
levels. The principal contribution to the change 
in the internal conversion coefficient comes, 
however, from the valence electrons, particu- 
larly the 4p, 4d, and 5s levels. In order to re- 
late compression to pressure, we have estimated 
the compressibility of technetium metal to be 
0.27 megabar~*. From this compressibility and 
the assumption that the internal conversion coef- 
ficient is linear in pressure, we calculate for 
the experimental pressure of 0.1 megabar a 
fractional decrease in lifetime of (2 to 4) x10, 
the quoted variation residing in the uncertainty 
of the structure of the 4p band. This result 
agrees with Bainbridge’s measurement, (2.30.5) 
x10°*, within the accuracy of our calculation. 


NUCLEAR LEVELS OF Cs***. M. G. Stewart and 
D. C. Lu, Institute for Atomic Research and De- 
partment of Physics, Iowa State University, Ames, 
Iowa (Received April 9, 1959; revised manu- 
script received December 28, 1959). 


The gamma rays of Cs*** following the electron- 
capture decay of Ba’** were studied using a coin- 
cidence scintillation spectrometer. Gamma rays 
with the following energies (in kev) were observed: 
386, 356, 301, 276, (220), 162, 82, 80, and 54. 
The energies (in kev) of the nuclear levels of Cs 
with their spins and parities are: ground state 
(7/2+), 82(5/2+), 162(3/2+), 383(3/2+), and 
438(1/2+). The electron-capture transitions and 
their relative intensities are as follows: (a) to 
the 438-kev level, 76%; (b) to the 383-kev level, 
11%; (c) to the 162-kev level, 13%. In the decay 
of Xe*’, (1.541.0)% of the beta transitions was 
observed to go to the 162-kev level in Cs*** and 
the remainder go to the 82-kev level. 


133 


INFLUENCE OF /] FORBIDDENNESS ON THE 82- 
kev TRANSITION IN Cs***, F. M. Clikeman and 
M. G. Stewart, Institute for Atomic Research 

and Department of Physics, Iowa State University, 
Ames, Iowa (Received April 9, 1959; revised 
manuscript received December 28, 1959). 


The gamma-gamma directional correlation of 
the 356-82 kev cascade in Cs*** was measured. 
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The directional correlation coefficients were 
A,/A, = 0.031 + 0.006 and A, /A, = -0.006 + 0.010, 
This is consistent with a spin assignment of 
3(Q)3(D + Q)3. The ratio of the E2 to M1 matrix } 
elements for the 82-kev transition is 6 = -0.139 
+0.007. The M1 transition is retarded by a fac. 
tor of ~700 and the £2 transition is enhanced by 
a factor of ~20, compared to single-particle 
transitions. The retardation of the M1 transitiq, 
is consistent with the assignment d,,. for the 82- 
kev state and g,,. for the ground state which wou 
make the M1 transition 7 forbidden. The enhance. 
ment of the £2 transition indicates there is a 
cooperative phenomenon present. 






TRITON ANGULAR DISTRIBUTIONS FROM 14.8. 
Mev DEUTERONS ON C*’. S. Mayo* and A. I. 
Hamburger, University of Pittsburgh, Pittsburgh, 
Pennsylvania (Received August 28, 1959). 


An elemental carbon target enriched to 66% in 
C!* was bombarded with 14.8-Mev deuterons. — 
Triton groups corresponding to the ground state, 
4.43-Mev, and 7.65-Mev levels in C” were mag- | 
netically analyzed and detected with photographic 
plates or a scintillation counter, at laboratory 
angles between 3 and 80 degrees. The angular 
distributions were found to be in good agreement 
with Butler pickup theory, for 7=1. The values 
of the cross sections at the peak of the angular 
distributions are 24.3, 9.6, and 0.36 mb/sterad, 
respectively. These values are accurate within 
+25%. The relative reduced widths connecting 
C** with the ground state, 4.43-Mev, and 7.65-Mev 
levels of C” are 1: 0.76: 0.039. The results are 


discussed. 
* 
Present address: Laboratorio Sincrociclotron, 
Comisi6n Nacional de Energia Atomica, Buenos Aires, 
Argentina. 
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RADIOACTIVE DECAY OF LUTETIUM- 174. 

R. G. Wilson and M. L. Pool, Department of 
Physics and Astronomy, Ohio State University, 
Columbus, Ohio (Received July 20, 1959; revised’ 
manuscript received September 17, 1959). 


Ytterbium oxide enriched to 98.4% in the 174 
mass number was irradiated with 6- Mev protons. 
An activity of approximately 165-day half-life 
was produced and assigned to Lu’™ by the identi- 
fication of the ytterbium K x-ray and of the ac- 
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tivities produced by similar proton irradiations 
of the other enriched isotopes of ytterbium. The 
observed activity of Lu’™ consists of the L and 
Kx rays of ytterbium and 76.6- and 1228-kev 
gamma rays which are in coincidence. Because 
no beta radiation exists in the activity of Lu’, 
the mode of decay is solely by electron capture 
to Yb'. Approximately 31% of the disintegrations 
of Lu’™* are to the ground state of Yb’™. In ad- 
dition to the 76.6-kev level of Yb'™, there is a 
1305-kev level with a spin of 0+. The transitions 
of Lu" to the 1305-kev level of Yb’™ are by L 
capture only and the percentages of electron cap- 
ture to the 76.6- and 1305-kev levels of Yb’ 

are approximately 59 to 10, respectively. A spin 
of 1- is assigned to the ground state of Lu’. 


CONNECTION BETWEEN NEUTRAL PION DE- 
CAY AND PROTON COMPTON SCATTERING. 
Maurice Jacob and Jon Mathews, Norman Bridge 
Laboratory of Physics, California Institute of 
Technology, Pasadena, California (Received 
August 20, 1959). 


The inclusion of the single-pion intermediate 
state significantly improves the agreement be- 
tween theory and experiment for proton Compton 
scattering. Excitation functions and angular 
distributions for various assumed values of the 
n° lifetime are presented. On this basis one may 
definitely exclude 7° mean lifetimes > 10~** sec 
or < 5x107!* sec. 


ACCURATE DETERMINATION OF THE i* 
MAGNETIC MOMENT. R. L. Garwin, * D. P. 
Hutchinson, S. Penman, and G. Shapiro, Co- 
lumbia University, New York, New York (Re- 
ceived August 4, 1959). 


Using a precession technique, the magnetic 
moment of the positive 1 meson is determined 
to an accuracy of 0.007%. Muons are brought 
to rest in a bromoform target situated in a ho- 
mogeneous magnetic field, oriented at right 
angles to the initial muon spin direction. The 
precession of the spin about the field direction, 
together with the asymmetric decay of the muon, 
produces a periodic time variation in the prob- 
ability distribution of electrons emitted in a fixed 
laboratory direction. The period of this variation 
is compared with that of a reference oscillator 
by means of phase measurements of the “beat 


note” between the two. The magnetic field at 
which the precession and reference frequencies 
coincide is measured with reference to a proton 
nuclear magnetic resonance magnetometer. The 
ratio of the muon precession frequency to that of . 
the proton in the same magnetic field is thus 
determined to be 3.1834+0.0002. Using a re- 
evaluated lower limit to the fnuon mass, this is 
shown to yield a lower limit on the muon g factor 
of 2(1.00122+0.00008), in agreement with the 
predictions of quantum electrodynamics. 


* 
Also at International Business Machines, Watson 
Laboratory, Columbia University, New York, New York. 


K* SCATTERING IN THE ENERGY RANGE 100- 
250 Mev. D. Fournet Davis, N. Kwak, and M. F. 
Kaplon, University of Rochester, Rochester, 
New York (Received August 21, 1959). 


Approximately 65 meters of K*- meson track 
were scanned in the energy interval 100-250 Mev. 
Analysis of the data shows that the K*-neutron 
charge exchange scattering rises appreciably 
above the value 1/5 corresponding to a pure 
J=1 interaction and at energies 2 150 Mev appears 
to be larger than the K -neutron direct scattering. 
Analysis of the angular distribution of the inelastic 
events for K* incident energies greater than 
150 Mev indicates a nonisotropic center-of-mass 
angular distribution, peaked towards 6, ., = 180°. 


RECOIL PROTON POLARIZATION FROM 225- 
Mev 1” -p SCATTERING. J. F. Kunze, T. A. 
Romanowski, J. Ashkin, and A. Burger,* Carnegie 
Institute of Technology, Pittsburgh, Pennsylvania 
(Received August 11, 1959). 


m~ mesons of energy 225 Mev were scattered 
from liquid hydrogen. The polarization of the 
recoil proton has been measured at two angles. 
For analyzing the polarization, a counter-con- 
trolled cloud chamber was used. Recoil protons, 
which experienced scatterings in a carbon plate 
of the cloud chamber, were photographed stereo- 
scopically. The pictures were later projected, 
and the necessary measurements of the scatter - 
ing were made directly in three dimensions. The 
computed polarizations are -0.13+0.16 ata 
laboratory recoil angle of 15°, and +0.36+0.29 
at 31°. The positive sign is for polarization in 
the direction of the vector cross product of the 
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incident pion momentum and the recoil proton 
momentum. A comparison of the data is made 
with various sets of scattering phase shifts which 
represent the differential cross-section data 
equally well. The results favor the Orear type 
of Fermi set in which the S-wave a, phase shift 
is positive. 

"Present address: CERN, Geneva, Switzerland. 


K~-MESON ELASTIC SCATTERING BY EMUL- 
SION NUCLEI. D. Fournet Davis, N. Kwak, and 
N. F. Kaplon, University of Rochester, Rochester, 
New York (Received August 14, 1959). 


Data on elastic nuclear scattering of K~ mesons 
in emulsion have been obtained. The differential 
cross section has been calculated by the WKB 
approximation in the method of partial waves. It 
is concluded that the real and imaginary poten- 
tials necessary to give the correct K~-nucleon 
cross section also give a good fit to the elastic 
nuclear scattering data. 


STOPPING OF 1~ AND K MESONS ON HYDRO- 
GEN IN NUCLEAR EMULSION. George A. Baker, 
Jr.,* University of California, Berkeley, Cali- 
fornia, and Los Alamos Scientific Laboratory, 
Los Alamos, New Mexico (Received August 26, 
1959). 


We show that the experimentally observed 
fact that the number of K” mesons stopping on 
hydrogen in nuclear emulsion is significantly 
larger than the number of 7~ is consistent with 
other known properties of these mesons. The 
reason is chiefly that the zero energy 1~ -p 
cross section is anamolously low. 


- 
Present address: Los Alamos Scientific Laboratory, 
Los Alamos, New Mexico. 


KINEMATICS OF GENERAL SCATTERING PRO- 
CESSES AND THE MANDELSTAM REPRESEN- 
TATION. T. W. B. Kibble,* Norman Bridge La- 
boratory of Physics, California Institute of 
Technology, Pasadena, California (Received 
April 27, 1959; revised manuscript received 
October 29, 1959). 


The kinematics of an arbitrary process in- 
volving two incoming and two outgoing particles 
is studied in terms of the invariants used in 
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Mandelstam’s representation, treating the three 
processes described by the same Green’s func- 
tion simultaneously. It is shown that the physica| 
regions for these processes are bounded by a 
cubic curve in the plane of the two independent 
invariants. The unitarity conditions are dis- 
cussed in the approximation of neglecting inter- 
mediate states of more than two particles. The 
formula for the spectral functions of the double 
dispersion relation is obtained explicitly in 
terms of the invariants chosen. 


*x 
Present address: Department of Mathematics, : 
Imperial College, London, England. { 


FOUR- FERMION INTERACTIONS WITH SPIN } 
NEUTRINO. C. L. Hammer and R. H. Good, 

Jr., Institute for Atomic Research and Depart- 
ment of Physics, Iowa State University, Ames, _ 
Iowa (Received August 19, 1959). : 


Beta decay involving a spin 3 neutrino in weak 
competition with the spin $ neutrino has been 
investigated. This process could explain the 


deviations from allowed shapes in the Gamow- f 


Teller interactions recently reported by Langer 
et al., but it would then be in disagreement with 
.-decay observations. 










INCONSISTENCY OF CUBIC BOSON-BOSON IN- 
TERACTIONS. Gordon Baym, Lyman Laboratory 
of Physics, Harvard University, Cambridge, 
Massachusetts (Received June 1, 1959). 


It is shown that there does not exist a ground 
state for a system of spin-zero Bose fields 
coupled only by local interactions involving three 
powers of the fields. Thus these interactions 
alone are not suitable for a model of interacting 
fields. 


REPRESENTATION OF STATES IN A FIELD 
THEORY WITH CANONICAL VARIABLES. 

F. Coester, Department of Physics, State Uni- 
versity of Iowa, Iowa City, Iowa, and R. Haag, 
Palmer Physical Laboratory, Princeton Uni- 
versity, Princeton, New Jersey (Received 
August 26, 1959). 


We investigate the properties of a functional 








representation of states for a self-coupled scalar 
field theory. The assumption is made that all 
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states can be generated by applying functionals 

of the field at a fixed time (¢=0) to the vacuum 
state. It is shown that for the class of models 
considered the Hamiltonian is uniquely deter - 
mined by the vacuum functional. The calculation 
of scalar products between states leads to func- 
tional integrals. The measure in this integration 
over function space is also determined by the 
vacuum state. Two methods for the evaluation 

of the functional integrals are discussed. The 
first one reduces the problem in some simple 
cases to the solution of an eigenvalue problem 

for a Hilbert-Schmidt kernel plus a finite number 
of ordinary integrations. The other one gives a 
perturbation series. 


SEVEN- DIMENSIONAL CHARGE SPACE. D. C. 
Peaslee, School of Physical Sciences, Austra- 
lian National University, Canberra, A. C. T., 
Australia (Received December 1, 1958; revised 
manuscript received October 12, 1959). 


Some implications of a charge space of seven 
(rotational) dimensions are considered for (I) 
particle classification, (II) strong interactions, 
and (III) weak interactions. The principal con- 
clusions from this point of view are that (I) all 
particles except the photon and graviton can be 
incorporated in 7-dimensional charge space; 
lepton conservation must be abandoned, which 
automatically introduces parity nonconservation 


into B decay; (II) the one-boson, two-fermion in- 
teraction is predominantly pseudovector in form 
and induces no mass differences; the 2-A mass 
difference arises from interference with a funda- 
mental two-boson, two-fermion interaction of 
lower charge symmetry; the =-N mass differ- 
ence has an “intrinsic” basis and is not due en- 
tirely to strong interactions; (III) the weak uni- 
versal Fermi interaction has signature -i under 
time reversal; the strangeness change AS =i 1 
is associated with 256 independent (in principle) 
terms; the rule AJ=4 + 1/2 needs further qualifi- 
cation because 1=T+ U, where T and U are in- 
dependent operators. 


THIRRING MODEL WITH VARIABLE INTER- 
ACTION. F. L. Scarf,* CERN, Geneva, Switzer - 
land (Received June 15, 1959). 


The Thirring model is solved with a variable 
coupling constant, A=A,f(x,¢). It is found that 
the infrared divergence is eliminated if f tends 
to zero along the past light cone. The S matrix 
is no longer diagonal in the physical particle 
representation and is generally not well defined 
in the sense of Haag’s theorem. The ordered, 
renormalized Heisenberg operators for », ~y* 
are computed and production processes are 
analyzed by examining matrix elements. 


a 
On leave from the University of Washington, Seattle, 
Washington. 











